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ABSTRACT 


The purpose of this thesis is to investigate the application of 
elastic theory to predict deflections in full-depth asphalt pavements. 
A section on the provincial Highway No. 15 northeast of Edmonton, 
Alberta, was selected for study. This pavement section was instrumented 
in two locations to record a) surface deflections and stresses induced 
at the pavement-subgrade interface, due to the application of standard 
18 kip single axle-loads, and b) various temperatures within the pavement 
structure. The resilient modulus of the subgrade soil was determined in 
the laboratory in a repeated loading triaxial test apparatus, whereas the 
stiffness of the asphalt concrete was obtained theoretically using nomo- 
graphic procedures. These properties of the soil and the asphalt con- 
crete were incorporated in the Chevron 5-layered elastic theory computer 
program and the surface deflections were computed in an iterative pro- 
cedure. The computed deflections were compared with the measured field 


deflections under similar environmental and test conditions. 


These comparisons showed that when the subgrade was divided into 
two layers, one for the compacted subgrade and one for the undisturbed 
subgrade, the computed and measured deflections were in reasonably close 
agreement. A design procedure was suggested utilizing the Chevron 
5-layered elastic theory computer program. This procedure takes into 
account the threeway interaction between tolerable deflection, perfor- 


mance and future traffic use as compiled by the Asphalt Institute. 
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CHAPTER I 


INTRODUCTION 


Pa eeoeneral 

There is a considerable interest on the part of highway engi- 
neers to develop a theoretical basis for designing and evaluating 
asphalt pavements. A number of methods under development involve the 
application of the theory of elasticity to layered systems as repre- 
sented by asphalt pavement construction. If it were possible to 
predict pertinent stress, strain and deflection responses due to wheel- 
loads by means of theory, a valuable step could be taken in the 
efficient design and evaluation of pavements and the economical utili- 


zation of pavement materials. 


Previous investigations in the application of the theory of 
elasticity to layered systems have been reported in the literature and 
are discussed in CHAPTER II. These indicate that elastic theory can 
reasonably be expected to predict deflections and strains due to the 
application of wheel-loads. This investigation is one of an increasing 
number of studies designed to verify the available solutions for pre- 


dicting strain and deflection in layered systems. 
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A major problem in solving various mathematical models for 
stresses and displacements in layered elastic systems by the use of 
computer programs, is in selecting the appropriate elastic moduli of the 
various pavement layers from subgrade to surface. For fine-grained 
cohesive soils commonly used as subgrade materials, the term resilient 
modulus has been proposed and is widely accepted for this purpose 


(Seed et al., 1967). 


Oo 
Mp os ao SOO Ceo OOO ce OOO Gao Oo Led 
x 
in which, 
Mp = modulus of resilient deformation in psi; 
Sq = repeatedly applied deviator stress in psi; and 
Set resilient axial strain corresponding to a 


particular number of stress repetitions in 


inches per inch (recoverable axial strain). 


1.3. Pavement Model 

A flexible pavement can be considered as consisting of layers 
of pavement material. Since the particular pavement under study in this 
investigation is a full-depth asphalt pavement, the materials contained 
therein are the asphalt concrete and the subgrade soil. The pavement 
model formulated for structural analysis consisted of a maximum of five 
layers so that an available elastic theory computer program could be 
utilized. Each layer can be visualized as consisting of nesta renes as 


material with very little variation in properties within it. 
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Various combinations in the number of asphalt concrete layers 
and the number of subgrade layers in the pavement model were tried, to 
find a model that would satisfactorily predict surface deflections and 
stresses induced on the subgrade due to wheel-loads when compared to 
those measured in the field. The stiffness modulus values assigned to 
the asphalt concrete were determined by the procedure suggested by 
McLeod (1969)* based on recovered properties of asphalt concrete and 
the resilient modulus values assigned to the subgrade soil were obtained 


from laboratory repeated loading tests on the subgrade soil. 


This is very similar to the approach used by Kasianchuk et al. 
(1969), except that they computed the tensile strain in the asphalt 
concrete to consider fatigue while surface deflections and vertical 


stresses in the subgrade were computed in this investigation. 


1.4. Objectives of the Thesis 

The primary objective of this thesis is to investigate the 
application of a particular elastic theory layered system analysis 
procedure to predict pavement responses and to compare these with field 
measurements. Surface deflections and vertical pressures measured at 
the pavement-subgrade interface were selected as the main responses for 


comparison. 


Since the "elastic" properties of the pavement layers are re- 
quired in the analysis, the resilient modulus approach was used to 


characterize the subgrade soil. The development of suitable laboratory 


* See List of References at end of text. 
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procedures for compacting and testing specimens at The University of 
Alberta has been a major phase of the investigation designed to meet 


this objective. 


l.og. Limitations of the Thesis 

The prototype pavement used in the study was a full-depth 
asphalt concrete pavement. Only two sections of this pavement were 
instrumented and tested during the fall of 1971; consequently data 
available for comparison were very limited. Hence, the findings of this 
investigation can only be applied to a full-depth asphalt concrete 


pavement. 


Of the many solutions available for predicting pavement re- 
sponses in layered systems, only the 5-layered elastic theory computer 


program developed by the Chevron Research Company was used. 


1.6. Organization of the Thesis 
The first chapter gives the introduction along with the pavement 
model, the objectives and the limitations of the thesis. The organi- 


zation of the thesis is also included. 


CHAPTER II reviews the literature on, a) Previous work, b) 
Repeated loading apparatus, which are important separate entities them- 


selves and must be taken into consideration in this investigation. 


The third chapter gives, a description Of the test projecu. 
Along with it is presented the data and test results obtained by the 
Research Council of Alberta, Highway and River Engineering Division. 
The data in this chapter are used in subsequent analyses while the 


deflection test results are used for comparison. 
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The fourth chapter outlines the testing procedure used, which 
constitutes a major part of the investigation. The exact operational 
procedure and the detailed description of various parts are given in 


APPENDIX A. 


The fifth chapter summarizes the test results and a detailed 


analysis is included. 


The sixth chapter gives the practical implications of the 


findings of this investigation and a design procedure is suggested. 


The seventh chapter is a summary of conclusions and recommen- 
dations drawn from the literature review and the results of tests and 


analyses. 


The Appendices follow a section on List of References and 
generally contain routine information on the subgrade material proper- 
ties, details of test results, details of apparatus and other infor- 
mation of minor importance to warrant its inclusion into the main body 


of the thesis. 
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CHAPTER II 
REVIEW OF LITERATURE 


2.0.) General 

This chapter reviews the literature necessary for this investi- 
gation and is divided into two major sections. The first section deals 
with previous work in the structural analysis of pavements and materials 
characterization of pavement components. The second section reviews the 


literature on various repeated triaxial loading devices. 


2.2. Previous Work 

The application of elastic theory to the design and evaluation 
of highway pavements has been widely favoured in the recent years. For 
more than two decades, a considerable amount of work has been done in 
this area of research most of which is reported in the Proceedings of 
the First and Second International Conferences on the Structural Design 
of Asphalt Pavements published in 1962 and 1967, respectively. A proper 
description of the behaviour of material under repeated loads, is indis- 
pensable for successful application of the elastic theory for layered 
systems to predict stresses, strains and deflections. The response of 
road building materials to repeated loading has been extensively studied 


by Seed and his co-workers at University of California, Berkeley, from 
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1955 and a summary of their findings has been reported by Seed et al. 


(1967). 


2.2.1. Approach to Structural Analysis 


Various multi-layer elastic theories have been proposed in the 
past, and these have been summarized by Whiffin and Lister (1962). The 
theories put forward are either two-layer or three-layer elastic 
theories. Some of these have been presented in tabular and graphical 
forms by Jones (1962), Peattie (1962), and Ahlvin and Uhlery (1962). 
With the advent of modern electronic computers, some of the computa- 
tional difficulties associated with the calculation of stresses, strains 
and deflections from the appropriate elastic theory for layered systems 
have been minimized. Thus, so long as appropriate material response 
characteristics are available, it would appear practicable to use such 
techniques to study the response of in-service pavements. The Chevron 
5-layered elastic theory computer program for calculating stresses, 
strains and deflections in a pavement has been successfully used by 
Monismith et al. (1967) and Shifley and Monismith (1968) in predicting 


pavement deflections. 


Various other researchers have demonstrated the applicability of 
elastic theory to predict stresses and deflections in a pavement. The 
theory has not only been applied to single-layer systems, but also has 
been extended to take into account multi-layered systems. A summary of 
the pertinent information from the previous studies is presented in this 


section. 
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Lister and Jones (1967) have examined in some detail the use of 
circular loaded areas to estimate the stresses existing in the different 
layers of the structural section. From their analyses they conclude 
that the assumption of a uniformly loaded circular area defined by the 
wheel-load divided by the inflation pressure will give results which in 
other than exceptional circumstances, are in error by less than + 2 per- 
cent when used to calculate interfacial pavement stresses and surface 
deflections under the whole range of tire-load conditions. They do 
note, however, that for high length to width ratios for tires which are 


associated with an overload condition, axial, radial and circumferential 


stresses will differ from those derived from circular assumption. 


Their analysis thus provides some confidence (at least to an 
engineering approximation) in the use of circular loaded areas, particu- 
larly since the Chevron 5-layered elastic theory computer program per- 
mits the determination of stresses and strains away from the axis of the 


load. 


Kirk (1967) has attempted to analyse the results of the AASHO 
road test within the framework of theory. For uncracked setedon®, he 
obtained a value of subgrade stress of one psi (based on fall deflection 
measurements) as being associated with the ability of the pavement 
sections to carry a large number of load repetitions for axle loads in 


the range,.12 to. 30 kips; 


Measurements of strains in actual pavements were reported by 
Gusfeldt and Dempwolff (1967), Klomp and Neisman (1967), and Nijboer and 


Delcour (1967). While the measured results and results computed by 
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9 
layered theory are not in precise agreement, all of these investigations 
demonstrate that layered system elastic theory may be quite useful from 
an engineering standpoint in assessing the response of pavements to 


various load conditions. 


Coffman (1967) has presented convincing data indicating that 
laboratory determined moduli for the various components of the pavement 
structure obtained from sinusoidal loading (complex modulus) can be used 
within the framework of layered system elastic theory to predict the 


surface deflection of pavements. 


Terrel and Krukar (1970) have used the Chevron 5-layered elastic 
theory computer program to predict deflections. Their predicted de- 
flections, using layered elastic theory along with material properties 
determined from dynamic tests on laboratory prepared samples, compared 
reasonably well with the measured deflections. They have also observed 
that although the values (like modulus of resilience) are accurately 
determined, it is difficult to assign appropriate values to these 
materials in the pavement at various times during the pavement's life. 
Hence, sampling and testing of the pavement materials obtained at the 
time of field measurements, would make it substantially easier to 
establish the effectiveness of this type of analysis and validity of 


elastic layer theory in general. 


Seed et al. (1967) attempted the correlation of surface deflec- 
tions predicted from repeated-load triaxial tests on representative 
samples with those measured in repeated plate load tests on a wide 


range of materials. They obtained reasonable agreement by assuming 
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10 
that the stresses were those predicted by Boussinesq theory and by 
taking a relevant soil modulus from triaxial tests carried out at con- 
stant radial pressure. Although both of these assumptions are definite 
over simplifications of the problem, the agreement observed by the 
authors suggest that within a limited range of test and boundary condi- 


tions, the errors are minor. 


They have also observed: 


"One of the greatest difficulties is presented by the 
variation of the resilient modulus of compacted clays with 
the magnitude of the deviator stress. Because the deviator 
stress induced in the subgrade due to any wheel-load appli- 
cation varies both vertically and horizontally, the subgrade 
in general will behave as a material with varying modulus, 
no matter how homogeneous it might be. Inclusion of such a 
variation in the elastic theory for layered systems has not 
been possible to the present time." 


From the preceding observation it is evident that the subgrade 
should be divided into several thin layers so as to have a minimal 


variation of resilient modulus within a layer. 


A number of solutions to boundary value problems of stresses and 
displacements in earth masses and layered systems are available and 
these are summarized in an annotated Bibliography (1969). Computer 
programs for these solutions are also available and Nair (1971) has 
pointed out: 


"Practicing engineers should note that computer programs are 
now available for solving layered system problems formulated 
in accordance with the methodology first outlined by 
Burmister. The number of layers these programs can handle 
covers the range of all practical problems. Because these 
programs have been 'debugged' and are operational, their use 
in routine design is a practical proposition. These programs 
provide information on stresses, strains and displacements 
throughout the pavement system." 
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One such computer program is a 5-layered elastic theory program 
developed by the Chevron Research Company, and this has been used in 


this investigation. 


2.202 Matorials Characterization 
a) Subgrade Soil 

Seed et al. (1967) have summarized their past investigations 
into the resilience characteristics of fine-grained cohesive soils 
(clays) under repeated loads. They have concluded that the factors that 
influence the resilience of clays under repeated loads are: a) Number 
of stress applications, b) Age at initial loading, c) Stress intensity, 
d) Method of compaction, e) Compaction density and water content, and 
f) Changes in water content and density after compaction. They have 
also demonstrated that kneading compaction produces laboratory specimens 
with resilience characteristics similar to those observed in field 


specimens for the same condition of test. 


b) Asphalt Concrete 


The stiffness of asphalt concrete must also be incorporated in 
the multi-layered elastic theory in order to investigate the response of 
a pavement to repeated loading. There are two ways of determining the 
stiffness of this and they are: 

1) Direct determination of mixture stiffness, in which 

stiffness is determined from repeated-load flexure or 


repeated-load compression tests. 


2) Indirect method of mixture stiffness determination, in 
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12 
which stiffness is determined from the recovered properties 


of asphalt concrete and a nomographic solution. 


The indirect method of mixture stiffness determination was used 
in this thesis, because of its quickness and simplicity. This method 
yields results that agree fairly closely with the corresponding results 
obtained from direct determination of mixture stiffness, as shown by 


Monismith (1966) and McLeod (1969). 


Indirect method of stiffness determination. Van der Poel (1954) 
was the first to suggest the term stiffness to represent the time and 
temperature dependent characteristics of asphalt concrete. This 


definition is represented by: 


S(t,T) = = of. agaress de ¢ open 
where, 
S = stiffness (in psi or kg/cm?) ; 
t = time of loading; 
T = temperature; 
o = axial stress (in psi or kg/cm?) ; and 
€ = axial strain. 


From a series of tests on both asphalts and mixtures he 
demonstrated that the stiffness of a mixture could be determined from a 
knowledge of the penetration and ring and ball softening point of the 
asphalt contained therein and the volume concentration of the aggregate. 
To simplify the stiffness computations, Van der Poel prepared 


nomographs. 
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Heukelom and Klomp (1964) have examined this approach in more 
detail and have suggested the following expression for determining the 
mixture stiffness using the same characterizing factors developed by 


Van der Poel. 


n 
Sani C 
. mix : Vee <2 3 - 2.2 
asphalt Vv 
where, 
de : B DN 
ae mixture stiffness (kg/cm*) ; 
= i 2 
Sasphalt stiffness of recovered asphalt (kg/cm) 
(from Heukelom and Klomp's nomograph, which 
is a slight modification of Van der Poel's 
original nomograph) ; 
op = volume concentration of aggregate defined as: 
volume of aggregate | 
volume of (aggregate + asphalt) ’ 
5 
n = 0.83 jg oe 
asphalt 


Heukelom and Klomp indicate that this expression is valid for 
mixtures with values for C,, ranging from 0.7 to 0.9 and air void 


CONEENts Ot the, order of S percent, that 1s, well compactéd mixtures. 


Van Draat and Sommer (1966) have suggested a correction to be 
applied to computed stiffness for mixtures with void contents in excess 
of 3 percent. This correction simply requires using a corrected volume 


concentration of aggregate in Equation 2.2. 


Cath wey OE UK ve See ee ts o's 200 


‘ 


m of dosovgas eids bonimexe sved (hORL) aqarolt bar moledyeh 


- 5 : s . niet) te — ih ie fn a. wi t Run 
(timvatsb ‘rot moreeotoxo gepkwollot sfy beteegaue over bag Tieseb 
ee 
VO y f E rid toarRHS eee ait3 gaiev 2e29mtite stwixnia 
wel wear 
-leod teb may 
( \ x, ' 
a 2 
: V Ge t i Kal ie 
. : a ee 
re : : ° 
: ; * * Mos 
J Ti sriqes 
\ i , " « ae : i rx i a , ‘ 
«A 
- s Le l my rw} te < o 
” 43 YO IST 2th e . _ , ) 
Li priqes 
, | OO 2 Bie mol IALIGI Cee | t)} 
wy 
j roi t H i r3 - a an ie | iff] c i> anf 
~ >fxigdersomon faentgrto 
Se 4 
at —_ 7) 
Pia 5 VG fF j iSsIMOD | Cav - Me 
5 : ¥ 
ra 
; ; 
e + 
Esti se . j (ae) Sirus fy 
COT yt : 
ea 7 gol cs. = nt 
i 
hay as 
SL RAG 2: 
ei a, . 
2 , wy © . : we in pofa bm ot - aay’ Os pope . 4 al . 7 
tot bilsv ef colezerqxe efit tad? otso thre Grell Das metexuen , 
iT ‘ y 
-< ‘ ni 
‘ } “ . ~ ry r re y 
DiGY 7 his 2,0 oF \/0 mort So ipagy .O noe esitav Ae tw 29 ar 
¢ 5 on 
; ae 
ba Ceeanktwene Te ; ee red ier - reeyt vi 
m@tTustkim bSjordtoO>S lilow ,.2eL JBAI .INISTEeQ. ZO Tepto sa zo. one J 
rh 
e — = ) 4 
, ne ; . t ——— ° AO TA. eran ; ee _ ~ - 
(ed of soitoerroo s betesague svad (ee!) temmee bas ta8%, as 


¥ oe 


% 


“ane ot bata .. 


te : as seri ” ‘ 
bert TH “Lomi ies nos itt smoot 
a ary Lory pais a2 


14 


in which, 
Cy = corrected volume concentration of aggregate; 
C, = volume concentration of aggregate defined in 
Equation 2.2; and 
AV = difference between actual air void content 


and the value of 3 percent (expressed in 
decimal form). 
It is suggested however, that the correction is applicable so 


long as the volume concentration of asphalt Cy exceeds $(1- Oe: 


McLeod (1969) has outlined a procedure for the indirect deter- 
mination of mixture stiffness based on recovered properties of asphalt 
concrete, and he found reasonable agreement with the measured dynamic 


modulus. 


Finn (1967) has explained the concept of the interchangeability 
of the time of loading and temperature for a mixture stiffness. He has 
also demonstrated how a single curve of stiffness versus time for a 
particular temperature (known as reference temperature) can be obtained 
for a particular mixture, knowing the stiffness versus time relationship 


for various temperatures. 


The methods suggested by Finn (1967) and McLeod (1969) have been 


used in this investigation. 


2.3. Repeated Loading Apparatus 


Highways carry mixed traffic at random intervals, moving at 
different speeds. Consequently, the pavement is rapidly loaded and 


unloaded when a vehicle passes over a particular section of it. This 
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15 
phenomenon is present during the pavement's entire life. Hence, it is 
desirable to consider repetitive loading when pavement materials are 


being evaluated. 


A repeated loading apparatus should ideally consist of a 
system that can load and unload samples at a desired frequency with 
negligible impact effect, and a system for measuring and recording the 
applied stress and resulting strains. Besides these requirements, the 
necessity for the measurement of associated parameters dictate other 


features of the apparatus. 


A number of different types of repeated loading apparatus have 
been developed by various investigators. In order to set up a repeated 
loading apparatus in the laboratory at The University of Alberta, it was 
considered necessary to review different features of the various 
apparatus for applying repeated loads, developed by other investigators. 
A summary of the essential features of other test equipment is pre- 


sented in TABLE II-l. 


The novel feature of the apparatus developed by Shackel (1970) 
is a hydromechanical device which can simulate a field stress path. A 
temperature controlled triaxial cell and a fully automatic data 


monitoring and acquisition system were incorporated in the equipment. 


The installation of a repeated loading triaxial test apparatus 
in the laboratory at The University of Alberta was under consideration 
for some time in the past. Progress had been made in the procurement 


of some parts thereby making the adoption of an apparatus modelled 
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after Seed and Fead (1959) feasible. Features of the apparatus are 


described in CHAPTER IV, with the details given in APPENDIX A. 


2.4. Summary 


The review of literature has disclosed that: 


1) 


2) 


3) 


4) 


5) 


The elastic theory can reasonably be expected to predict 
stresses and deflections due to the application of wheel 


loads. 


The 5-layered elastic theory computer program developed 
by the Chevron Research Company is one of the successful 
computer programs available for solving practical problems 


in layered systems. 


The resilient modulus is found to satisfactorily represent 
the elastic modulus of the subgrade soil for use in these 


computations. 


The stiffness of asphalt concrete mixtures can be 
determined theoretically by the indirect method, according 


to the procedure suggested by McLeod (1969). 


An apparatus, such as that developed by Seed and Fead 
(1959), is suitable for repeated-load triaxial tests on 


subgrade materials. 
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CHAPTER III 


TES? "PROJECT 


Jala) 16sec) Prolect and ite 

The development of an elastic design method for pavements has 
been under investigation for some time under the Cooperative Highway 
Research Program, comprising the Highway and River Engineering Division 
of the Research Council of Alberta, the Department of Highways and 
Transport of Alberta, the Water Resources Division of the Department of 
the Environment of Alberta, and the Department of Civil Engineering at 
The University of Alberta. An extensive study of the elastic properties 
of a section of full-depth asphalt concrete pavement on a clay subgrade 
at several locations on a four-lane divided freeway section has been 
undertaken by the Highways Division of the Research Council of Alberta, 
with this particular laboratory investigation as a part of the 


Cooperative Program. 


The test project was the new Highway 15-A-1, Alberta, which is 
located northeast of Edmonton. It starts from the city limits at mile 
7.11 and ends just east of the junction of Highway No. 37 at mile 
15.19. This highway is a full-depth asphalt concrete pavement con- 
structed on a clay subgrade. The highway subgrade was constructed in 


the late summer of 1971 and the test sections were completed in the 
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19 
fall of the same year. The location of the instrumented sections on 


this highway are at mile 9.44 and mile 11.00 as indicated in FIGURE 3.1. 


The asphalt concrete was laid directly on the surface of the 
compacted subgrade, in three layers with nominal thicknesses of 3, 3 and 
2 inches, respectively. The total thickness of the asphalt concrete 
pavement according to the Contract and Specifications of the Department 
of Highways and Transport was 8 inches, but due to construction varia- 
tions it was 7.2 inches at mile 9.44 and 9.3 inches at mile 11.00. The 
cross-sections of the highway at these two locations are shown in 


PIGHRES S.2 and 3.5, respectively, 


The asphalt used as the asphalt concrete pavement binder was 
specified as an AC 275 with a minimum viscosity of 275 poises at 140°F 


and a penetration of 250 at 77°F. 


3.2. Instrumentation 

During construction of the highway thermocouples and pressure 
cells were installed within the pavement structure by the staff of the 
Research Council of Alberta, Highways Division. These installations 
were at mile 9.44 and mile 11.00, in the inner lane of the westbound 


lanes on the outer wheelpath. 


Standard Copper-Constantan thermocouples were installed at these 
locations at approximately 2 inch depth increments from the pavement 
surface within the asphalt concrete pavement, and thereafter at 2 feet 
intervals within the subgrade to a depth of 96 inches below the pavement 


surface. 
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The pressure cells installed in the pavement were fabricated in 
the laboratory, by the staff of the Research Council of Alberta, 
Highways Division. This was a half-bridge system with an active gauge 
in one arm of the bridge and a compensating gauge in the other arm of 
the bridge. The active gauge was a spiral strain gauge while the 
compensating gauge was a quarter inch metal foil gauge. These gauges 
were encased in an aluminum, cylindrical hollow box, 4.5 inches in 
diameter and 0.5 inches thick, so that the active gauge was placed 
against the top plate, which also acted as a diaphragm. The pressure 
cell was connected to a Sanborn 321 recorder. At one particular 
location, three such pressure cells were embedded in the subgrade such 
that the tops were flush with the subgrade surface and were 24 inches 


apart along the outer wheelpath. 


3.3. Measurements 

After the second layer of asphalt concrete was laid, static 
moduli of the pavement and the subgrade were determined from the 
Benkleman Beam Test according to the procedure described by Shields 
and Hutchinson (1961). The dynamic moduli of the pavement and the 
subgrade were determined by the Vibrator Technique according to the 
procedure described by Shields (1971). Vertical stress on the subgrade 
under the wheel-load as registered by pressure cells, the temperatures 
as registered by thermocouples, and surface deflections were measured 
using the Benkleman Beam Test according to CGRA (1965) rebound pro- 
cedure... The» tests, wexearepeated-atter itheatinalmlife of asphalt 


concrete was laid. Surface deflections recorded are shown in FIGURE 


# 


ti & Sitdst stow tremevay oAY af boliptent efles saueeotq of 
st2 ofa yd .vrodexodel odd 
itos we Activ moteye sybivd-tisd s eeweitl! .noieivad eyewigen 


) Mts TAATO ot Ee oaE ifjsamoq@nod s bik ogbiad od7 to. mts eno 


\ " * ‘ e. ee 
_ 1 ca lie Pa, | SD niee et oe Pee seur a6 : PO em — pr ere Fy 
xa MG wail : PEPE rhe | t (BISN HOMl LST TAtp & 248W sues Bas ~ SAn@amo af 


7 


4 
r PR hy ral 1m Seber PP ae pede +4 lin ‘| hase o 
irc. & ,xod wollon [BorYbRelvo ris As AI bezahons sie, 


vi i if Vt t t J oct xX pfs pao | biit G 0) EE | & 
. 4 4 = on P nto ies oov ads 
i } oA l ee t s th Dai ’ A J yaNY J Iti 
fi ;  ~SbToo9% Ike RECOURSE B * perp eliiels 
: 7 . a rf 
53 OSphbs Was, 2 a a rq fine SeTny , 
TOW. DE e shavvaue ent ctwenenit sraw 206 
4 P 
isclesiiw ratyo srt 
f 
j 
esas motuenoM - 
ee en ny 
f enw 33 v ) SiBATSS 160 Tavs baocasée soit pexra ae 


ns moxt poninxeteh srew sbsredue ons Bris o_o toh 
DLS £1 ‘d beditaesh exubesotqg oly of yaubrooga Jeer nso8 na 
/ 


oy t thoneysy sat to fighom simengb on - (het) no neaiediat bas 


Sat Ot gribrasos ep indosT toferdi¥ eat vo boretm teh ozow db 


Sbhstgdue os no ezsvt2 lsotitev -€IVOL) abl sine pean ie 


—=— 


* 

ea Sa egmes oz cpa “A 

Lae a % ie : 
bh basa ee 

e 


ane a 


24 
3.4 a), b), c), and d) as longitudinal surface deflections, while the 


measured vertical stresses and moduli are tabulated in TABLE III-1l. 


Tests on recovered asphalt concrete were performed to evaluate 
the properties of the mixture and the asphalt used. Penetration, 
viscosity, and ring and ball softening point values of the asphalt; 
and dry density, bulk specific gravity, percent air voids, percent 
volume of mineral aggregate (VMA), and the percentage of asphalt content 


in the mixture were determined and are tabulated in TABLE III-2. 


These measurements were used subsequently for analyses and 


comparison with theoretically predicted values. 


The pressure cells were calibrated by embedding them in sand 
and recording the output which was then compared with calculated 
stresses by means of Boussinesq theory. As such the estimated accuracy 
of the pressure cells are + 15 percent. This is sufficient to give an 
idea of the magnitude of stresses induced in the pavement. The 


accuracy of the thermocouples are + 1.5°F and that of the dial gauges 


in the Benkleman Beam are + 0.001 inches. 
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TABLE III-1 
FIELD MEASUREMENTS TAKEN FROM HIGHWAY 15-A-1 


Pavement Modulus (psi) Pavement Modulus (psi) 
Temperature Temperature 


fe he bce Fecal | 
9.44 27.2 60.0 75.0 269300 14600 57.0 61.0 1410000 
WBEL,OWP,IL 11.8 42.0 42.0 454000 13600 42.0 42.0 1700000 
11.00 i 18.4 42.0 44.0 300000 14000 
WBL,OWP,IL 11.9 36.0 39.0 464400 14200 


Pavement 
Thickness 


Vertical 
Stress 


subgrade 


psi 


TABLE III-2 - 
RECOVERED PROPERTIES OF ASPHALT CONCRETE FROM HIGHWAY 15-A-1 


Location Dry “| Bulk Asphalt Air Voids Penetration at} Viscosity Ring & Ball 
Density Specific Content at 60 C Softening 
Gravity Point 


4 


NOTE : VMA : Voids in Mineral Aggregate 
C, : Volume Concentration of Aggregete 
cy : Corrected Volume Concentration of Aggregate 
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CHAPTER IV 


TEST EQUIPMENT, PROCEDURE AND PROGRAM 


4.1. General 

This chapter summarizes the testing program and the testing 
procedures followed in this investigation. A description of the re- 
peated loading triaxial test apparatus and associated equipment is 
given. Methods of sample preparation and compaction of samples are 


briefly described. 


4.2. Repeated Loading Triaxial Test Apparatus 


The primary objective of the testing program was the determi- 
nation of the resilient modulus of the fine-grained cohesive subgrade 
material under various deviator stresses. A repeated loading triaxial 
test apparatus was needed for this purpose and the equipment designed 
and constructed after the model of Seed and Fead (1959) was put into 


operation as part of this investigation. 


A steel framework made of channel sections was built for 
seating three triaxial cells along with their loading devices and 
allied equipment required for their operation. All the three loading 
devices were designed to produce the same load characteristics, so that 


consistent results could be anticipated. 
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Load is applied to a specimen in a triaxial cell by means of 
a piston and loading yoke, the load being repeatedly applied to the yoke 
by means of an air-pressure cylinder mounted below the triaxial cell. 
For this purpose compressed air is passed through a regulator and stored 
at a constant pressure in a reservoir tank. The tank is connected 
through a solenoid valve to the pressure cylinder. When the valve is 
open, air pressure is transmitted through the cylinder and loading yoke 
to the specimen. When the valve is closed, the air in the pressure 
cylinder passes through an exhaust pipe in the valve and the load is 


removed. 


The load applied to the specimen is varied by regulating the 
air pressure in the cylinder. By calibrating the load against the air 
pressure as recorded by a pressure gauge, the desired load can readily 


be obtained. 


Application of pressure to the air-pressure cylinder is con- 
trolled by a solenoid-operated three-way valve. Movement of the valve 
to admit air or exhaust air from the pressure cylinder is controlled by 


an electrical timing unit. 


The specimen in the triaxial cell is surrounded with water. 
Air pressure from the reservoir tank is applied through a pressure 
supply valve in the base of the triaxial cell, which constituted the 
confining pressure and is transmitted to the specimen through the con- 
fining water. The confining pressure is recorded by a pressure trans- 
ducer, which is connected to a strain indicator. Knowing the 
calibration factor of the pressure transducer, the strain indicator 


reading is converted to obtain pressure. 
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The specimen is loaded axially by placing the entire cell under 
the loading yoke. Deformation of the specimen is measured by a LVDT, 
which is clamped to the piston of the triaxial cell applying the load 
and the moving end is screwed to the cell cap. The magnitude of the 
load is measured by a load cell on the loading yoke which sits with a 


firm contact on the piston of the triaxial cell. 


The LVDT and the load cell are connected to an ultra-violet 
recorder through a galvanometer and signal conditioner unit. The loads 
and deformations are recorded on a light-sensitive paper in the ultra- 
w2olet recorders) Ihe calabration)ofsthe LVDT and the load cell is so 
matched with the galvanometer and signal conditioner unit, that the 
magnitude of deformation and load is read directly from the recording 


paper. 


The number of load applications to a specimen is recorded 
automatically in a six-digit resettable counter connected to the 


control valve. 


The triaxial cell used is a standard size suitable for testing 
cylindrical samples 4 inches in diameter by 8 inches high, as described 
by Bishop and Henkel (1962), except that the recess under the base was 
adjusted so that the cell could sit on the seat in the frame. A 
threaded hole on the cap of the triaxial cell was made for the LVDT to 


be screwed to the cap. 


PLATES 1 and 2 show the test apparatus with further details 


given in APPENDIX A. 
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PLATE «1. 


A) GENERAL VIEW 


RECORDING EQUIPMENT AND TIMING UNIT 


REPEATED LOADING TRIAXIAL APPARATUS AND 
EQUIPMENT USED IN TESTING PROGRAM 
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CLOSE-UP OF AN ASSEMBLED TRIAXIAL CELL 
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Aes. Preparation of Soil Samples 


The size and shape of the samples to be prepared, as determined 
from the test apparatus, were essentially 4 inches in diameter and 8 
inches high cylindrical specimens. Specifications of the Department of 
Highways and Transport of Alberta required the pavement subgrade to be 
compacted to one hundred percent Standard Proctor Density at optimum 
moisture content. Hence it was decided to maintain the same requirement 


for compaction of laboratory samples. 


The embankment material of the prototype pavement was 
obtained from borrow pits in the field. This material was spread on a 
flat surface and allowed to air-dry. After air-drying it was crushed 
and passed through a No. 4 sieve. Material passing through this sieve 


was used for sample preparation. 


The Standard ASTM D698-70 Test for determining moisture-density 
relation of soils, was performed with the soil and the optimum moisture 
content and the corresponding dry density was determined. Hygroscopic 
water content of the soil was determined and distilled water was added 
accordingly to obtain the desired moisture content. One half percent 
extra water was added to allow for evaporation during mixing. Soil 
was mixed thoroughly by hand for five minutes. Hand mixing of the soil 
was done, because of the large amount (4,000 grams) of soil required for 
the preparation of one sample, which could not be accommodated in the 
mechanical mixer available in the laboratory at The University of 


Alberta. Material for six samples was mixed at the same moisture 
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35 
content and placed in a plastic bag. The mouth of the bag was tied 
tightly and the bag was stored in the moisture room for twenty-four 


hours to allow for moisture migration within the soil. 


4.4. Compaction of Samples 


auAn) Impact Compaction 


At the beginning of the investigation, the kneading compactor 
was not available. As a result, samples were hand compacted using a 
drop-hammer hand compactor. The compactor was a 5.5 pound drop-hammer 
with a 12 inch free fall, as specified for Standard ASTM D698-70 Test. 
The mould used was a 4 inch diameter and 8 inch high split mould with 
a 2 inch collar at the top. Trial compactions were performed to 
determine which procedure could be used to attain the Standard ASTM 
D698-70 Densities. It was found that soil compacted in five layers 
with twenty-five blows per layer was close to the required density. To 
control the uniformity of samples, material required for preparation of 
one sample was divided into five equal parts. Each part was placed in 
the mould and spread evenly by hand after which twenty-five blows of 
the drop-hammer were applied around the mould. After the last blow to 
the fifth layer, the collar was removed and the excess material was 


trimmed. 


4.4.2. Kneading Compaction 


The kneading compactor was installed in the laboratory at The 
University of Alberta and was made available for use at a very late 


stage of the investigation, hence only a limited amount of work could 
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34 
be done with samples prepared by kneading compaction. The kneading 
compactor was calibrated by trial and error so as to attain Standard 


ASTM D698-70 Densities for the particular sample size. 


From the experience gained it was observed that Standard ASTM 
D698-70 Densities were attained by compacting the sample in five 
layers, with twenty-five tamps per layer and a ram dwell time of 1.5 
seconds. For soils at optimum moisture content, this density was 
achieved by applying a foot pressure of 150 psi. However, when the soil 
was wet of optimum moisture content, this foot pressure was too high and 
caused excessive penetration of the foot (about one and one half inches). 
Consequently, the foot pressure was reduced to 100 psi, which gave a 
foot penetration of about one quarter inch. The details of calibration 


are included in APPENDIX B. 


PLATE 3 shows the kneading compactor with further details 


given in APPENDIX B. 


4.5. Outline of Testing Program 

The testing program consisted of repeated loading triaxial tests 
on compacted clay samples to determine the resilient modulus (Mp) under 
various deviator stresses. It was decided to use clay samples at 
optimum moisture content and wet of optimum moisture content to see how 
the resilient modulus would be affected under adverse conditions of 
construction. It was also decided to see how different methods of 
compaction affected the resilient modulus. For this purpose, three 
series of tests were run. One series of samples were prepared with 


impact compaction, one series with kneading compaction, and one series 
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with field cores obtained with Shelby Tubes. The results of the tests 
were presented as points on curves of Mp versus deviator stress. In 

some cases a single specimen was tested more than once, an acceptable 
procedure according to Terrel (1972), provided that the stress levels 


and the number of repetitions were not excessive. 


All the samples were subjected to repeated loading at the rate 
of 20 load applications per minute, with a load duration of 0.25 
seconds. The testing program followed in this investigation is 
summarized in TABLE IV-1. The program is divided into three main 
series: a) tests on impact compaction samples (L series), b) tests on 
subgrade cores (F series), and c) tests on kneading compaction samples 


(K series). 


Available literature indicated that a confining pressure of 
3 psi was usually used in determining the resilient modulus of subgrad 
soil in the laboratory. This magnitude of lateral pressure is usually 
exhibited in the subgrade near the pavement-subgrade interface, and is 
dependent on the thickness of the pavement. Thus a confining pressure 


of 3 psi was adopted for use in this testing program. 
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CHAPTER V 
RESULTS AND DISCUSSION 


Se aeeenera | 

This chapter is divided into three major sections. In the first 
section, curves are presented from which the stiffness modulus of 
asphalt concrete can be determined theoretically. The second section 
consists of the subgrade materials characterization, a part of which is 
derived from the results of the laboratory repeated loading triaxial 
tests on the subgrade soil. The third section uses the moduli values 
provided by the first two sections and involves structural analyses to 
find an effective pavement model for predicting surface deflections and 
stresses induced on the subgrade, with the help of the Chevron 5-layered 


elastic theory computer program. 


hia Asphalt Concrete Stiffness 


The top layer of a flexible pavement usually consists of 
asphalt concrete whose stiffness varies with temperature and loading 
time as noted by Van der Poel (1954). Such being the case for the 
problem under consideration, it was imperative to determine the stiff- 
ness of the asphalt mixture under conditions of varying temperatures 


and loading times. 
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The stiffness of the asphalt mixture used in the pavement was 
determined by an indirect method based on the recovered properties of 
the asphalt concrete (TABLE III-2) and nomographs presented by 
Heukelom and Klomp (1964), using the procedure suggested by McLeod 
(1969). The indirect method and its development is discussed in 
CHAPTER II. A master curve to obtain stiffness of the asphalt concrete 
at various loading times for a reference temperature of 60°F was deter- 
mined and is presented in FIGURE 5.1. In order to use this master 
curve over a range of temperatures, the shift factor concept, as de- 
scribed by Finn (1967) was used and a variation of shift factor with 
temperature is shown in FIGURE 5.2. Sample calculations for the 
preparation of FIGURES 5.1 and 5.2, along with a typical example of 


using these figures are included in APPENDIX D. 


Finn (1967) has shown that for a material exhibiting thermo- 
rheologically simple behaviour, the position of the modulus curve but 
not its shape, is altered on the time scale by temperature changes. 
That is, if one were to define the modulus for a wide range of times at 
a series of temperatures, these curves at jdifferent: temperatures, are 
shifted only along the time axis. If a reference temperature T, is 
chosen then the amount by which another temperature curve is shifted 
along the time axis so as to superimpose on the reference temperature 


curve, is known as the shift factor and is physically defined as: 
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42 
where, 
Ar = eshitt ‘factor; 


ty = time required to observe a phenomenon at 


temperature T; and 


t, = time required to observe the same phenomenon 


at some reference temperature T). 


This procedure helps to reduce the family of curves at different 
temperatures to a single curve for a particular reference temperature. 
In addition, the plot of shift factor versus temperature enables the 
determination of shift factor for any intermediate temperature within 
the range of plotted temperature, thereby reducing many additional 


computations. 


5.3. Subgrade Soil 
The subgrade soil under the pavement was an inorganic clay with 
medium plasticity and can be classified as CL type of soil in the 


Unified Classification System. The characteristics of the soil are: 


Liquid Limit = 38.5% 
Plastic Limit = 16.4% 
Shrinkage Limit = 14.3% 


ASTM D698-70 TEST: Maximum Dry Density 2.2 YHs7 £te 


Optimum Moisture Content = 15.2%. 


Samples prepared in the laboratory were tested immediately 
after their preparation. The first set of three samples of the L test 
series were subjected to a minimum of 25,000 load repetitions, detailed 


calculations of which are included in APPENDIX C. It was observed from 
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43 
the results of this series that the resilient modulus did not change 
Significantly after 5,000 load applications. Thus it would be 
reasonable to assume that the resilient modulus at 5,000 load repeti- 
tions is representative of the material under the particular test 
conditions. Although other investigators have chosen to compare results 
at 25,000 load repetitions, this shorter test procedure appears justi- 


fied, in view of the preliminary testing. 


The results obtained from the laboratory testing program of the 
subgrade soil are summarized in TABLE V-1. This table gives the resili- 
ent moduli at 5,000 load repetitions, of samples subjected to various 
deviator stresses, along with their water contents and dry densities. 
The estimated errors in the resilient modulus determination varied from 
0.9 to 14.3 percent, the magnitude of error at a low deviator stress 
being high. However, the majority of errors were in the range of one to 
four percent. A sample calculation of the error is given in APPENDIX A. 
On a close examination of the table it will be seen that as the deviator 
stress increases the resilient modulus decreases and so also does the 


error. 


The subgrade soil tested in the laboratory under repeated loading 
was divided into three test series designated as L, F, and K test series. 
Series L comprises samples compacted by impact compaction, the F test 
series were Shelby Tube cores obtained from the field, and the K test 
series were samples compacted by kneading compactor. The results of the 
three test series given in TABLE V-1 are plotted separately to show the 


variation of resilient modulus with stress intensity. 
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TABLE V-1 
SUMMARY OF LABORATORY REPEATED-LOAD TEST RESULTS 


Sample Water Dry Deviator Contining Resilient Modulus a¢ Escimcted 
Number Content Density Stress Pressure 5000 Load Repetitions Error in M. 


pcf psi psi 
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Note : Series L - Impact Compaction Standard ASTM Density - 112.2 Ibs/ft? 
Series F - Shelby Tube Cores OMC - 15,2 % 
Series K - Kneading Compaction 


* Suspected malfunction of the apparatus; repeated test with L 7,8,9,10,11,12. 
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5.3.1. Laboratory Prepared Samples 


FIGURE 5.3 shows the variation of resilient modulus with 
deviator stress for the L test series. The results were comparable to 
those observed by Monismith et al. (1967) and Seed et al. (1967) for a 
Similar material. It is seen that for samples near the optimum moisture 
content, the rate of decrease in the resilient modulus was very pro- 
nounced with a slight increase in the deviator stress especially in the 
low range of deviator stresses (up to 8 psi), while in the higher range 
of deviator stresses, the rate of decrease in the resilient modulus with 
increase in the deviator stress was low. It is also observed that for 
samples wet of optimum moisture content, the pronounced decrease in the 
resilient modulus in the low deviator stress range is not present and 
the rate of decrease in the resilient modulus with increase in the 
deviator stress over the whole range of deviator stresses is moderate. 
In addition to this, the curve as a whole shifts to a lower order of 
resilient modulus, indicating that increasing the moisture content 
above the optimum results in a high resilient deformation which is 


undesirable in a pavement structure. 


The samples of the L test series were prepared by hand com- 
paction and preparation of identical samples having similar densities 
was difficult. Hence, the results obtained from testing these samples 
would be expected to show some variation. It is seen that for this 
type of compaction, a variation of one percent in moisture content and 
three pounds per cubic foot in dry density within a set of samples, does 


not affect the resilient modulus significantly. 
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The effect of stress intensity on resilient modulus for K test 
series is shown in FIGURE 5.4. The trend of the curves are similar to 
those presented in FIGURE 5.3 and discussed earlier. The samples were 
mechanically compacted resulting in more uniform and closer to identical 
samples being obtained (TABLE V-1). Except for points (1) and (2), the 
results fit in a smooth curve with low scatter. Another test was run 
with the same test conditions as was used for samples which yielded 
points (1) and (2). The new results obtained were plotted in FIGURE 
5.4 and it was found that points (1) and (2) shifted to points (1A) and 
(2A), respectively, to fall on the curve indicating that there was some 


error in the first test. 


A foot pressure of 75 psi was used while preparing sample K4. 
The foot penetration was about one quarter inch and the dry density 
achieved was lower than that of the Standard ASTM D698-70 Test for the 
same moisture content. Consequently, the foot pressure was increased 
to 100 psi for preparing samples K5 and kK6. The foot penetration in 
this case was about one half inch and the dry density was increased to 


that of the Standard Density at that particular moisture content. 


It would be expected that at identical moisture contents, 
samples with lower densities would yield lower resilient moduli than 
samples with higher densities. But FIGURE 5.4 indicates that sample 
K4, with a lower density than those of KS and K6, gave higher peeuiaent 
moduli than those for K5 and K6 in the same range of deviator stresses. 
It is rather interesting to note that the amount of foot penetration 


during compaction of samples effects the resilience characteristics of 
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49 
clayey soils and the results here show that samples prepared with a 
lower foot penetration gave a higher resilient modulus than those pre- 
pared with a higher foot penetration. The effect of foot penetration on 
the resilience characteristics of soils, also noted by others warrants 


further investigation but is considered beyond the scope of this thesis. 


There was no significant effect on the resilient modulus when 
samples were tested more than once. It can be seen from FIGURES 5.3 and 
5.4 that the test results fit reasonably well in the curvilinear rela- 
tionship with no appreciable scatter. Since a single sample was used up 
to four times in these tests, the results show that there is little 
effect on the resilient modulus by repeated use of the same sample, 
which was also observed by Terrel (1972). However, this aspect was not 
verified by a detailed investigation, but from the results of this 
investigation there is reason to believe that a sample may be used 
several times to produce compatible results, provided that the applied 
stress, the resulting strain and the number of load repetitions are 


kept low. 


Seoee. ~6Field Samples 


FIGURE 5.5 shows the variation of resilient modulus with 
deviator stress for the F test series. Only one set (comprised of 
three samples) was used in the testing program and the results of these 
tests show a scatter of points and no definite curve can be drawn 
through these points. However, most of these points could well be in- 
cluded in a band, with a curve being drawn through the upper and lower 


series of points. If this band is treated as representative of the 
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on 
variation of the resilient modulus with the deviator stress for this 
test series, it can be seen that the trend of this band is similar in 
nature to that of the curves in FIGURE 5.4, with the upper and lower 
boundaries resembling the 15.1 percent and 18.2 percent water content 


curves, respectively. 


The samples of the F test series as stated earlier in this 
chapter, were Shelby Tube cores and these were subjected to some 
disturbance while sampling. The recovered length of the samples were 
less than 8 inches (between 6 and 8 inches), making it difficult to 
maintain the length to diameter ratio of two, as in the case of labo- 
ratory prepared samples. These could possibly be one of the reasons 
responsible for the.wide scatter of test results in this series, but it 
is considered that the primary reason for the scatter resides in the 
fact that of the three samples tested, no two samples possessed 


identical moisture contents and dry densities, as revealed in TABLE V-1. 


From the experience thus gained by testing the field samples, 
it is evident that great care must be exercised while obtaining field 
samples. In future, if block samples are obtained rather than cores, 
the size of the block being large enough so that at least three samples 
could be carved from it, the samples would be subjected to a minimal 
disturbance while sampling. Therefore, it could be expected that test 
samples carved from the block would possess identical properties and 


the test results would be more consistent. 
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5.4. Structural Analysis 

A difficult part in predicting stresses and deflections in a 
pavement system using the multi-layered elastic theory computer pro- 
grams, is the selection and assignment of proper material properties of 
the various layers from subgrade to surface of the pavement. Terrel 
and Krukar (1970) found it difficult to assign appropriate moduli values 
to the pavement materials at various times during the pavement's life 
and have suggested to use the pavement properties at the time of field 
measurements, which would then make this type of analysis effective and 


establish the validity of the elastic layer theory. 


The pavement under study was a full-depth asphalt concrete 
pavement and the materials contained therein are the asphalt concrete 
and the subgrade soil, whose proper moduli values are essential for 
the effective use of the elastic theory. The modulus of the asphalt 
concrete was obtained from the curves presented in FIGURES 5.1 and 5.2. 
Laboratory resilient modulus test results (FIGURES 5.3, 5.4, and 5.5) on 
the embankment material were accepted as representative moduli values 
for the subgrade soil. Pavement models were formulated which contained 
a maximum of five layers in order to be used with the Chevron 5-layered 
elastic theory computer program. Various combinations in the number of 
subgrade and asphalt concrete layers were used to arrive at a pavement 


model which would yield reasonable predictions of pavement responses. 


Since the predicted pavement responses were compared with the 
pavement responses obtained from the CGRA Benkleman Beam Test, a wheel- 


load of 9,000 pounds with a tire pressure of 80 psi was used in the 
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53 
structural analysis. Although the actual loading was with a vehicle 
having an 18,000 pound axle loading with dual tire configuration, the 
simplifying assumption of single rather than dual contact areas has 
been made. From available literature an appropriate value of Poisson's 
ratio for asphalt concrete was assumed to be 0.33 and that for the sub- 
grade soil between 0.4 and 0.5. The modulus values were assigned to the 
various layers and the Chevron 5-layered elastic theory computer program 


was used to compute the stresses and defléctions. 


5.4.1. Single-Layer Pavement 

The pavement system was analysed to see the effect on stresses 
and deflections by dividing the subgrade into layers. The results are 
presented in TABLE V-2 and the deflection profile for each case is 
shown in FIGURE 5.6. During field measurements at mile 9.44 it was 
observed that the temperature of asphalt concrete in the pavement was 
uniform from top to bottom at 42°F. In the analysis, the asphalt 
concrete was maintained as a single layer and the stiffness assigned to 
it from FIGURES 5.1 and 5.2 was that which corresponded to a temperature 
of 42°F and a loading time of 0.1 seconds. The subgrade was treated as 
a single layer with semi-infinite thickness and the effect of assigning 


low and high moduli values were investigated. 


In one case, a modulus value of 10,000 psi was assigned to the 
subgrade. This modulus value corresponds to a deviator stress of 
approximately 12 psi (FIGURE 5.4), which is likely to be expected in 
the subgrade at the pavement-subgrade interface for this pavement 


thickness. The computed center point stresses and deflections are 
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TABLE Y-2 


MILE 9.44, PAVEMENT TEMPERATURE 42°F, 


Layer Thickness Poisson's Average Layer Mo: 


Temperature 
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LOADING TIME 0.1 
SECONDS 


duli and Stiffness of Pavement ( psi ) 
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( psi ) 
Measured parameters under the wheel : a) Surface Deflection : 20.9 x rs ae ins. 
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55 
shown in column (A) of TABLE V-2. In another case, a modulus value of 
25,000 psi was assigned to the subgrade. This modulus value corre- 
sponds to a deviator stress of approximately 1.5 psi (FIGURE 5.4), which 
is likely to be expected in the subgrade at a depth of three feet. The 
computed center point stresses and deflections are shown in column (B) 


of TABLE V-2. 


It is observed from the results presented in columns (A) and (B) 

of TABLE V-2 and the corresponding deflection profiles in FIGURE 5.6 
that there is a great disparity in the stresses and deflections for 

the two cases. However, center point: stresses and deflections in 
column (A) are closer to the measured parameters indicating that when 

a two-layered system is used for the structural analysis, the modulus 
value assigned to the subgrade must be in the region corresponding to 
the likely deviator stress in the subgrade at the pavement-subgrade 


interface. 


In subsequent cases, the subgrade was divided into two, three 
and four layers with the last layer always being of semi-infinite 
thickness and the preceding layers adding up to a thickness of 66 
inches, which is the thickness of compacted subgrade in the prototype 
pavement. The moduli values assigned to the various layers were taken 
from the experimental results (FIGURE 5.4) corresponding to the esti- 
mated deviator stresses in those layers. The center point stresses 
and deflections computed for these cases (C, D, and E) were in very 
close agreement with the measured values, the maximum variations being: 
a) 1.5 x 1073 inches in deflection, which is about 7 percent, and b) 


0.5 psi in stress, which is about 4 percent. 
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Examination of FIGURE 5.6 shows that deflection profiles corre- 
sponding to the cases in columns (A), (C), (D), and (E) of TABLE V-2, 
are Similar in nature to the measured profile and they follow it very 
closely. The deflection profile for case (C), which is a three-layered 
system, has the closest agreement with the measured profile, the 
difference of deflection under the wheel being two and one half percent. 
The dividing of the subgrade into further layers causes the surface 
deflection to decrease, but the rate of decrease is moderate. In 
comparing with the measured deflection under the wheel, the deflection 
for a five-layered system (case E) is underestimated by seven percent 
while for a two-layered system (case A) it is overestimated by twelve 
percent, although all the computed deflection profiles are almost 


parallel to each other. 


The results presented in TABLE V-2 and FIGURE 5.6 and their 
close agreement with the respective measured parameters is encouraging, 
and seems to justify the use of elastic layer theory for multi-layered 
pavement systems. However, it must be noted that the conclusion thus 
drawn is based on the utilization of pavement properties at the time of 


field measurement, and an assumed time of loading. 


5.4.2. Two-Layer Pavement 

Another analysis similar to the previous case was done. The 
asphalt concrete was divided into two layers instead of a single layer 
as in the previous case, and a reasonable temperature difference of 15°F 
between the top and bottom of the pavement was assumed. After several 


trials (TABLE V-7) to determine the appropriate loading time which will 
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57 
be discussed later, a loading time of 0.5 seconds and the average 
temperature of each layer were used to determine the stiffness of each 
layer, while the modulus assigned to the subgrade was determined as in 
the previous case. The subgrade was divided to constitute a maximum 


of three layers. 


The results of this analysis are shown in TABLE V-3 and FIGURE 
5.7. The stresses and deflections obtained are almost identical with 
those obtained in the previous analysis. Since the subgrade moduli in 
both cases are in the same range, it is the stiffness of the asphalt 
concrete which influences the magnitude of stresses and deflections. 
Hence proper use of the stiffness values are of great importance for 
predicting stresses and deflections. If the asphalt concrete within 
the pavement has a uniform temperature, it can be treated as a single 
layer and the assignment of stiffness values is fairly simple. On 
the other hand if there is temperature variation within the asphalt 
concrete, attempts must be made to obtain the actual pavement stiff- 
nesses. The behaviour of asphalt concrete is complex and its stiffness 
varies with time and temperature as noted by Van der Poel (1954). Under 
these circumstances obtaining actual stiffnesses is difficult and great 


care must be exercised in choosing the effecting parameters. 


In this analysis it was found that a combination of loading 
time of 0.5 seconds and temperatures of 38°F and 30°F for the top and 
the bottom layer, respectively, is equivalent in terms of stiffness to 
a combination of a loading time of 0.1 seconds and a uniform temperature 


of 42°F throughout both the layers. These stiffnesses can be said to be 
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COMPUTED STRESSES AND DEFLECTIONS FOR A TWO-LAYER PAVEMENT 


MILE 9.44, SURFACE TEMPERATURE 40°F, PAVEMENT THICKNESS 7.2 INS., LOADING TIME 0.5 


SECONDS 


Layer Thickness Poisson's Average Layer Moduli and Stiffness of Pavement ( psi ) 


Temperature 


175000 
380000 


SURFACE DEFLECTION 
UNDER THE WHEEL 


( x 1073 


VERTICAL STRESS ON 
SUBGRADE UNDER THE 
WHEEL ( psi ) 


Measured parameters under the wheel : a) Surface Deflection 
b) Vertical Stress on Subgrade : 11.8 psi. 


175000 175000 
380000 380000 
10000 10000 
15000 20000 


25000 25000 


: 20.9 x 107° ins. 


WHEEL LOAD POSITIONS (INS.) 


0 10 20 30 


10 


15 


SURFACE DEFLECTION (x107? INS.) 


20 


NOTE : 


25 


--" 
- 
-°* 


ToS secneos, 


so 


--7 
- 
== 


MEASURED DEFLECTIONS 
A : THREE LAYER SYSTEM 
B ; THREE LAYER SYSTEM 
C : FOUR LAYER SYSTEM 
D : FIVE LAYER SYSTEM 
E :; FIVE LAYER SYSTEM 


THIS GRAPH IS TO BE READ IN 
CONJUNCTION WITH TABLE V-3 


FIGURE 5.7. COMPUTED LONGITUDINAL DEFLECTION PROFILES 
FOR A TWO-LAYER PAVEMENT 
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59 
appropriate representations in the pavement model since in both cases 
the stresses and deflections computed were in very close agreement with 


the measured values. 


5.4.3. Variable Pavement Thickness 

It has been seen earlier that division of the subgrade into two 
layers, one layer for the compacted subgrade and one layer for the 
undisturbed subgrade, produced results which were the closest in 
agreement with the measured parameters. Hence, structural analyses with 
a Similar model were made for different asphalt concrete thicknesses for 
which measured stresses and deflections were available for comparison. 
The results of the analyses are presented in TABLES V-4, V-5 and V-6. 
It.is seen that when_the asphalt. concrete is treated as a single layer 
as in the case for 5.5 inches pavement thickness, the stresses and 
deflections computed and shown in TABLE V-4 compare favourably with the 
measured values when a loading time of 0.1 seconds is used. When the 
asphalt concrete is treated as two layers as in the case for 9.3 inches 
pavement thickness, the computed parameters shown in TABLE V-5 compare 
favourably with the measured parameters when a loading time of 0.5 
seconds is used. This phenomenon of 0.1 seconds for single-layer 
pavements and 0.5 seconds for two-layer pavements was also observed in 


the two previous analyses. 


Computed stresses and deflections for a 5.1 inches thick 
pavement given in TABLE V-6 are not in precise agreement with the 
measured values. It is observed from the previous analyses that the 


temperature of asphalt concrete was low (in the region of 40°F) and the 
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TABLE V-4 
COMPUTED STRESSES AND DEFLECTIONS FOR 
5.5 INCHES THICK PAVEMENT 


MILE: 11.00, SURFACE TEMPERATURE: 45°F, PAVEMENT THICKNESS: 5.5 INCHES 


Layer Thickness Poisson's Average Layer Moduli and Stiffness at Loading Times (Secs.) of: 
Ratio Temperature (in psi) 


5.5 . 250000 185000 155000 135000 121000 
12.0 . 10000 100090 10000 10000 10000 
Semi Infinite . 25000 25000 25000 25000 25000 


SURFACE DEFLECTION 
UNDER THE WHEEL 


(x 1073 ins.) 


VERTICAL STRESS ON 
SUBGRADE UNDER THE 
WHEEL ( psi ) 


Measured parameters under the wheel : a) Surface Deflection 222.5, xlOm ins. 


b) Vertical Stress on Subgrade : 18.4 psi. 


TABLE V-5 
COMPUTED STRESSES AND DEFLECTIONS FOR 
9.3 INCHES THICK PAVEMENT 


MILE: 11.00, SURFACE TEMPERATURE: 40°F, PAVEMENT THICKNESS: 9.3 INCHES 


Layer Thickness Poisson's Average Layer Moduli and Stiffness at Loading Times (Secs.) of: 
Ratio Temperature (in psi) 


360000 260000 217000 190000 175000 
380900 280000 230000 205600 185000 


10000 10000 10000 10000 
25000 25000 25000 25000 


SURFACE DEFLECTION 
UNDER THE WHEEL 


(x 107° ins.) 


VERTICAL STRESS ON 
SUBGRADE UNDER THE 
WHEEL ( psi ) 


Measured parameters under the wheel : a) Surface Deflection : 16.5 x eee ins. 


b) Vertical Stress on Subgrade : 11.9 pai. 
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62 
stiffnesses used in the analyses gave good predictions. But in this 
case, the temperature within the pavement varied from 75°F at the top to 
60°F at the bottom, which could be a possible reason for the predicted 


values to be inconsistent with the measured values. 


The asphalt used in the prototype pavement was of a low vis- 
cosity and high penetration grade, with a penetration index of -1.0 
indicating that this asphalt is highly temperature susceptible. This 
type of asphalt is known to exhibit a complex behaviour at high tempera- 
tures (Kopvillem and Heukelom, 1969), thereby resulting in a possible 
erroneous determination of the mixture stiffness at 75°F. As discussed 
im section 5.2,mthe concept of shift factoreas observed by Finn (1967) 
is only applicable to materials exhibiting thermorheologically simple 
behaviour. As such, the shift factor used for determining the stiffness 
at higher temperature was physically possible from FIGURE 5.2, but the 
actual effect was not known. When using the concept of shift factor, 
one should be aware of its possible limitations and the range of its 


applicability. 


5.445) | Variacionsin Loading Time 


The effect of variation of loading time from 0.1 to 10.0 seconds 
on the computed stresses and deflections was investigated and the re- 
sults are tabulated in TABLE V-7. The subgrade was treated as two 
layers and the asphalt pavement as two layers. The subgrade moduli were 
kept constant and the pavement stiffness varied with loading time. It 
was observed from the results that with this increase in loading time, 


the computed stresses increased approximately three-fold and the 
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64 
deflections increased approximately two and one half times over the 
whole range under consideration, the rate of increase for each 


increment of loading time being moderate. 


It was also observed from this analysis that when the pavement 
was treated as two layers with a temperature variation within it, a 0.1 
seconds loading time no longer predicted precise pavement responses as 
in the case for a single-layer pavement. These trials indicate that a 
loading time of 0.5 seconds with the combinations of temperatures 
considered in each pavement layer gives an appropriate stiffness, which 
when used in the pavement model predicts reasonable stresses and 
deflections. The trials in this analysis are the basis of selecting 


the loading time for the analysis in section 5.4.2. 


5.4.5. Variation in Temperature 

Initially, a hypothetical pavement section was chosen to formu- 
late a pavement model to see the effects of variations in a) pavement 
temperature from 40°F to 70°F, and b) the thickness and moduli of 
compacted subgrade, on computed stresses and deflections. The results 
of these analyses are tabulated in TABLE V-8. The pavement model 
consisted of three asphalt concrete layers having a total temperature 
variation of 15°F between the surface and the bottom of the asphalt 
concrete pavement. The subgrade was divided into two layers, one for 
the compacted subgrade and one for the undisturbed subgrade. The 
loading time was kept constant at 0.1 seconds and the stiffness varied 
with temperature. At a particular temperature, the pavement stiffness 


was kept constant and the compacted subgrade modulus was varied while 
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66 
the modulus of the undisturbed subgrade was maintained at 25,000 psi. 
The first trial incorporated a compacted subgrade thickness of 66 
inches and the second trial incorporated a compacted thickness of 24 


inches. 


It was observed from the results of this analysis that as the 
pavement temperature increased, thereby decreasing the pavement 
stiffness, the computed stresses and deflections also increased. At 
a particular temperature and subgrade modulus, the decrease in sub- 
grade thickness resulted in a decrease in deflection and an increase in 
stress. Also, at the same temperature and subgrade thickness, an 
increase in the subgrade modulus resulted in a decrease in deflection 
and an increase in stress. This is what would be naturally expected 
because an increase in modulus means a stiffer material and hence it 


is likely to undergo less deformation. 


It can be concluded from this analysis that at higher tempera- 
tures (70°F), the thickness and moduli of the compacted subgrade have 
considerable effect on the computed stresses and deflections, while at 
low temperatures (40°F) the effects are moderate. The increase in 
temperature increases the computed pavement responses, the stresses 
increasing approximately three-fold and the deflections increasing 
approximately two and one half times over the whole range of tempera- 
ture (40°F to 70°F), the rate of increase between the increment of 60°F 


and 70°F being quite pronounced. 
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5.5. Summary 


The major observations in this chapter can be summarized as: 


1) 


2) 


3) 


4) 


5) 


6) 


7) 


The resilient modulus of a clayey soil does not change 
Significantly after 5,000 stress applications, when the 


samples are tested immediately after their preparation. 


The rate of decrease in the resilient modulus was very 
pronounced in the low deviator stress range, and it was 


low in the higher range of deviator stress. 


The resilient modulus versus deviator stress curve shifts 
to a lower order of resilient modulus when the moisture 


content in the soil is increased beyond optimum. 


The foot penetration of the kneading compactor affects 


the resiizence of clayey soils. 


A sample can be tested several times, under repeated 
loading, to produce compatible results provided the 
applied stress, the resulting strain and the number of 


load repetitions are low. 


When proper material properties are incorporated in the 
Chevron 5-layered elastic theory computer progran, 
reasonably good prediction of pavement responses are 


obtained. 


When the subgrade is divided into two layers, one repre- 


senting the compacted subgrade and one representing the 
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undisturbed subgrade soil, and the modulus for each layer 
is assigned from the laboratory tests, best correlation 
between the predicted and measured deflections are 


obtained. 


8) The behaviour of asphalt concrete is very sensitive to the 
temperature and loading time and care must be taken in 


considering these factors when stiffnesses are determined. 


The subgrade moduli obtained from the field static test using 
the deflection bowl analysis and presented in TABLE III-1l, are within 
the range of the resilient moduli determined in the laboratory for the 
same material. This provides some confidence in using the deflection 
bowl analysis to arrive at the subgrade moduli, and either method could 


be used to yield reasonable results. 


oe 
P i F ’ 


, “+ ’ ry * \ = : 
Nel fone TOP @ulubow oft bans »1 foe abprgdy2 bodrutseibmy 


yf heed er shelinn co “cbse aoe ana b lt Sweet ee be 
iif s&s 35, D> Tes ,&s2o3 vrata todel one mers i MIZLe@2ea 2k 


— Ds Naas hal aves ek dele ia, oT Pe | ais mel 
OTH “MOL TISLISD bowers BNR SSIDLoOSIC ony AeowsoG 
P ’ : 
, —— 
re) 
j . Pe." : ~ a 
manrevTrTar 
» Wo IZBTCS <a 
on’ 
> pun yey opt ay i otces » es le 7 ye 7 
} OVIPLeN@e Ye 8797 2HOD tiedqzse 0 twolyeded ef (@ 


. 
5 sll ¢ 4 _ ~ 6, dl ae oe ~ 7." 4 Vored = 
L it i Jeum oreo Dig omit gaibeol bre onijoreqess 


- 


, ¢ - ‘ am +" es ; oat se ws ae ere rime? P = 
mMoIao i 28e2ONI77ite Wow eyoFos? gears SMT LINLEaNOS 
oe. 
4 ce. a a oe a ee par i: sf" Fl 7 
[J t OAS MOT? SORLSEIGO tlLiubom S58 radue ent . 


iu 
{ 


ry J 
a 


tiw ete ,! | SIUAD af bataseetq Bis ebevyiens Iwod mottze{2eb | 


— : 7 
ach ‘elnatt es Oe ae ee ae Ye ‘a } (kai ele Sx a 
af3 Lot yroteTousi suit nt benimveteb Elubom daeklifteeyt eAy to Sunset s 
bs 42 . 
; ’ 
aI it 3 2nebiines emoe eobivorg afAT .f 
4 £ . : ne a ebaees ; oe ee eee oe » > gla ei ie ag eee 
DLO bonyom fo bas ,tigbom BhatTgdve si3 1g evigis of eteyisna 


% 


t sifsnoeset bisty ot beew 
* 


CHAPTER VI 
PRACTICAL IMPLICATIONS 


Gr. General 

The previous chapter indicated that when proper pavement 
materials properties were used in the 5-layered elastic theory, reason- 
ably good predictions of pavement deflections and stresses were obtained. 
The asphalt concrete stiffness was very sensitive to temperature and 
loading time. This latter parameter is difficult to assess, however, if 
several full-depth sections for which data are available for comparison © 
are analysed at a particular temperature, the loading time can be 


deduced with reasonable accuracy. 


The elastic theory for multi-layered pavement systems may be 
useful in developing design curves for full-depth asphalt pavements. A 
detailed investigation in this aspect was not undertaken, but a pro- 
cedure illustrating the development of a method and possible use is 


outlined. 


6.2. —lolerable Defléctiens 

Many studies in the past have been made to relate pavement 
deflection with performance and future traffic use, the major findings 
of which have been summarized by Kingham (1969). He has compiled a 


traffic-deflection relationship shown in FIGURE 6.1, from the experience 
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71 
of various agencies in different countries. The Asphalt Institute has 
adopted a Design Line which is considered to represent a conservative 
relationship in that the probability of an unsatisfactory design is very 
low. Detailed investigation into this aspect of the problem was not 
undertaken because it is considered to be a study by itself. However, 
the consideration of the Asphalt Institute's Design Line for use with 
the design curves suggested here was checked with the data available in 


this investigation. 


6.3. Design Curves 

It has been indicated in the previous section that a guide such 
as the Design Line of the Asphalt Institute is available to relate 
pavement deflection with performance and future traffic use. This 
information can be effectively used for design purposes if a relation 
between the pavement deflection and pavement thickness can be estab- 
lished. The Chevron 5-layered elastic theory computer program which 
computed pavement deflections that agreed reasonably well with the 
measured values in this investigation, can be used for this purpose. 
A procedure for developing design curves using this computer program is 


outlined in this section. 


For simplicity in developing the design curves, a two-layered 
pavement model with one layer of asphalt concrete and one layer of sub- 
grade, is used here for illustration. If more accuracy is desired two 
layers of subgrade (one for compacted subgrade and one for undisturbed 
subgrade) can be used in the pavement model to be incorporated in the 


computer program. 
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Initially, for a particular subgrade modulus and asphalt 
concrete stiffness, the Chevron 5-layered elastic theory computer pro- 
gram is used to compute the surface deflections of the pavement for 
asphalt concrete thicknesses in increments of two inches. Without 
altering the subgrade modulus, asphalt concrete stiffness is changed to 
the next value and pavement deflections are computed for the range of 
asphalt concrete thicknesses used in the previous case. This process 
is repeated for the entire range of pertinent asphalt concrete 


stiffnesses. 


The surface deflections of the pavement computed for different 
asphalt concrete thicknesses and stiffnesses at that particular subgrade 
modulus, are plotted in a relationship of surface deflection versus 
pavement thickness. This will result in curves for a particular 
stiffness extending over the range of asphalt concrete thicknesses used. 
A typical plot of this nature for a subgrade modulus of 10,000 psi is 
shown in FIGURE 6.2. Similar curves are obtained for each different 


subgrade modulus in the range of the appropriate subgrade moduli. 


6.4. Use of Design Curves 

One of the current criterion for the design of pavements and 
overlays is to keep the peak pavement deflections within certain 
tolerable limits (Kingham, 1969). It has been seen in section 6.2 that 
maximum tolerable deflections of a pavement can be selected from the 
experience of various agencies in different countries, by knowing the 
volume of traffic the pavement section is expected to carry. The sub- 


grade is the weakest in spring and peak deflections of pavements are 
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74 
most likely to be obtained at that time; hence, the spring environmental 
conditions should be used in choosing the pavement materials properties 


and responses. 


The maximum tolerable deflection the pavement is allowed to 
undergo is selected from FIGURE 6.1, using the Asphalt Institute's 
Design Line and the expected volume of traffic. The subgrade modulus 
for spring conditions is obtained and a set of design curves similar 
to the one shown in FIGURE 62? but corresponding to the determined 
subgrade modulus is chosen. The value of the tolerable deflection is 
entered in the deflection axis and is followed horizontally to the 
estimated spring pavement stiffness curve (interpolation may be re- 
quired) from where the corresponding thickness of asphalt concrete 


pavement is obtained. 


This technique of determining asphalt concrete thickness may 
also be used for the design of overlays of existing pavement structures 
by formulating a suitable pavement model to be used in the computer 
program. The pavement model thus formulated can be checked by compari- 
son with a few measured Benkleman Beam deflections and corrections to 
the model if necessary can be implemented. For a desired deflection, 
the difference between the thickness thus obtained and the existing 


thickness of asphalt concrete is the thickness of overlay. 


The consideration of the Asphalt Institute's Design Line for 
use with the design curves suggested here was checked with the data 
available in this investigation. It was found that when working back 


from the design curve (FIGURE 6.2) using the materials properties in the 
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fall and converting the fall deflection to the peak deflection in 
spring (by multiplying the fall deflection by 1.8, as observed by 
Kingham), the Design Traffic Number arrived at was approximately sixty- 
five. This, according to the Asphalt Institute's (1969) traffic 
classification represents medium traffic, which is reasonable for the 
highway under study in this investigation. Hence, the Asphalt 
Institute's Design Line is reasonable for consideration in obtaining 


tolerable deflections for use in the design method suggested here. 


Considering practical aspects such as variations during con- 
struction and the capability of construction equipment, it would seem 
reasonable that the thickness be rounded off to the next half inch. 
If such is to be the case, a two-layered pavement model is considered 


sufficiently accurate to estimate the required thickness. 
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CHAPTER VII 
CONCLUSIONS AND RECOMMENDATIONS 


@-13 'CONGLUSTONS 

The main purpose of this thesis was to investigate the effec- 
tiveness of the elastic theory for layered systems in predicting 
pavement responses. The conclusions reached from this investigation 
were drawn from the repeated loading triaxial test program, the analyses 
of data using the Chevron 5-layered computer program and the literature 
review. The conclusions presented are limited by the scope of the 
testing program and availability of data for comparison. The pavement 
under study was a full-depth asphalt concrete pavement, hence the 
findings of this investigation can be only applied to that type of 


pavement. 


a) The conclusions drawn from the review of literature are as 
follows: 

1) A proper description of the behaviour of pavement materials 
under repeated loads is indispensible for successful appli- 
cation of the elastic theory for layered systems. For this 
purpose, the resilient modulus is found satisfactorily to 
represent the elastic modulus of the pavement materials for 


use in computations. 
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The 


the 


1) 


2) 


3) 


Vins 


Kneading compaction produces laboratory specimens with 
resilient characteristics similar to those observed in the 


field specimens. 


Asphalt concrete stiffness can be determined by the use of 
nomographic procedures based on the recovered properties of 


asphalt concrete. 


conclusions made from the results of the testing program on 


subgrade soil are: 


The determination of the resilient modulus at 5,000 stress 
applications when samples are tested immediately after 
their preparation in the laboratory, and using a sample 
several times under different test conditions provided the 
applied stress, the resulting strain and the number of load 
applications are kept low, has been found to be a satis- 


factory procedure in this investigation. 


The effect of stress intensity on the resilient modulus of 
clays is such that the rate of decrease in the resilient 
modulus is very pronounced in the low deviator stress range 
and it is almost constant in the higher range of deviator 
stress; and the curve as a whole shifts to a lower order of 
resilient modulus when the moisture content is increased 


beyond the optimum. 


The foot penetration of the kneading compactor affects the 


resilience of clays. 
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78 
It is concluded from the analyses that, 


1) When the subgrade is divided into two layers, one repre- 
senting the compacted subgrade and one representing the 
undisturbed subgrade, and the modulus for each layer is 
asSigned from laboratory tests, the Chevron 5-layered 
elastic theory computer program computes stresses and 
deflections which are in very close agreement with the 


measured values. 


2) The behaviour of asphalt concrete is very sensitive to the 
temperature and loading time and care must be taken in 


considering these factors when stiffnesses are determined. 


3) From the construction point of view, a two-layer pavement 
model incorporated in the Chevron 5-layered computer 
program is sufficiently accurate for the design of full- 


depth asphalt pavements and overlays. 


Recommendations 


The apparatus used in this investigation performed satisfactorily 
and gave excellent results. However, the electrical timing unit 
sometimes failed to operate during the test and its replacement 
is desirable. The timing unit has another disadvantage in that 
it can give a minimum on-load time of 0.25 seconds, which 
represents the loading time of a vehicle travelling at approxi- 
mately 15 miles per hour. When the timing unit is being re- 


placed, it would be advantageous to have some extra features 
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79 
included in the new one. It is recommended that the new timing 


unit have the capability of, 


1) Achieving the on-load times of vehicles travelling at 


high speed. 
2) Simulating the loading sequence of multiple axle vehicles. 


3) Operating continuously during the test for any number of 


load applications. 


It was found in this investigation and from available literature 
that the foot penetration of a kneading compactor affects the 
resilient modulus of clayey soils. It was also observed from 
the literature review that kneading compaction produced 
laboratory samples with similar resilience characteristics as 
those in field samples. In view of the importance of kneading 
compaction, a detailed investigation into the effect of foot 
penetration on the resilience of different types of soils would 


prove useful for future works. 


Most highways are using treated bases for their construction. 
An investigation similar in nature to that in this thesis, could 


bescarriecd OuLeto, ditferent types of treated bases: 


The procedure for obtaining design curves arrived at in this 
thesis, should be checked with more of the existing full-depth 


asphalt pavements to strengthen its validity. 
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APPENDIX A 


REPEATED LOADING TRIAXIAL TEST APPARATUS 


Awl. Description 


The repeated loading triaxial test apparatus set up at the 
University of Alberta is modelled after the one developed by Seed and 
Fead in California. FIGURE A.1 shows the layout of the apparatus. The 
essential features and the functions of the various components are de- 
Scribed in this” section. “FIGURE A.I° should be referred to for corre- 


sponding alphabetic letters. 


A Counterbalance for Loading Yoke. As the name implies, it is 
a weight used to counterbalance the effect of the weight of the loading 


yoke on the piston and ram. 


B Load Cell. A device for measuring the deviator stress 
applied to the specimen. This replaces the conventional proving ring, 
because of its ability to be used with the electrical recording system. 


Operating range : 0 to 2,000 lbs. 


C LVDT. Linear variable differential transformer, commonly 
known as LVDT, is a displacement transducer operated electrically for 
measuring the deformation of a sample. Here it is used to measure the 
axial deformation of the sample external to the cell D. The LVDTs used 
here are Hewlett Packard displacement transducers, model 7DCDT-1000. 


Displacement range : + 1.000 inch (full scale). 
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D Straaxtalecell.” This triaxial*cel? 1s of a standard size for 
testing 4 inch diameter and 8 inch high cylindrical samples. The base 
is supported by a platform in the testing machine frame, only over its 
central area so that the loading yoke can pass by the base without 
obstruction. The removable cylinder is of a 0.375 inch thick perspex. 
Because of the butt joint in the perspex, which is the source of 
weakness, four circumferential bands of one inch width and 0.1 inch 
thick and made of resin-bonded glass fiber, are fixed on the outside of 
the cylinder as reinforcements. The loading ram is a 0.75 inch diameter 
ground stainless-steel rod, which is fitted with a collar and a hemi- 
spherical piece. The features of the triaxial cell are shown later in 
this section. Capacity : Loading ram - 6,000 lbs. axial load, 


Cell - 150 psi. 


E Bay Operation Switches. These switches when switched on, 


start the operation of the respective bays. 


F Guide Frame for Loading Yoke. This unit consists of a plate 
with a slot bolted to the frame, and an annular cylinder threaded in the 
inside and smooth outside. The cylinder is screwed to one arm of the 
loading yoke and is made to pass through the slot of the plate. Such a 
unit is installed on each arm of the loading yoke. The purpose of this 
unit is to restrict the movement of the loading yoke to a vertical 
direction thereby insuring that the load applied to the load cell is 
vertical. This unit thus helps in preventing the load cell from being 


damaged. 
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G Loading Yoke. A rectangular frame, connected to the piston 
at the bottom and to the loading ram of the triaxial cell through the 
load cell at the top. It transfers the load from the piston to the 


loading ram. 


H Pressure Transducer with Three-Way Valve Switch. The pres- 
sure transducer is a very sensitive gauge. It is used here to measure 
the confining pressure in three triaxial cells. (A single transducer is 
used to measure the confining pressures of three cells, individually.) 
This is achieved by using a four-way valve in which the pressure trans- 
ducer is connected to one of the arms while the other three arms are 
connected to the three triaxial cells. By using a three-way valve 
switch, the pressure from two triaxial cells is shut off and the 
pressure from the third triaxial cell is measured by the transducer. 


Capacity : 300 psi (range, 0 to 300 psi) for transducer No. 28749. 


J .Air-Pressure Cylinder and Piston. This unit transfers the 
air pressure from the air-pressure reservoir tank to the loading yoke, 
in a mechanical form. The ones used here are : Hannifin Series 2A, 
Style HB (NFPA Style MF 6) with 4 inch diameter bore and LB+Stroke = 


4.875 inches; Rod and Thread - KK, Style 4 and 9 (0.75 inch - 16). 


K Air-Pressure Regulators. These units are installed at 
various tappings on the air-pressure reservoir tank for supply of air 
pressure to the triaxial. cells. These units regulate the air pressure 


to the triaxial cells. Capacity : 125 psi maximum operating pressure. 
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L Safety Microswitches. These switches are installed on the 
frame close to one arm of the loading yoke, for each bay. A steel bar 
is sandwiched between two nuts and is mounted on the arm of the loading 
yoke nearest to the microswitch. The position of the steel bar on the 
arm of the loading yoke can be adjusted by rotating the nut on which 
the steel bar is mounted, thereby adjusting the gap between the switch 


and the steel bar. 


The maximum axial deformation of the sample is estimated and the 
position of the steel bar is adjusted accordingly. In case there is a 
failure of the sample, the steel bar will come down and hit the micro- 
switch. This will cut off further operation of that bay. Hence, damage 


to the apparatus is avoided. 


M Counters )This*is® aysix-digit’ resettable electrical: counter, 
which is connected to the solenoid valve, and records the number of load 


applications on the specimen. 


N Air-Pressure Reservoir Tank. This is a cylindrical tank 
which stores compressed air under pressure. The pressure requirements 
for the three triaxial cells can be supplied by this single tank. 


Capacity : 125 psi maximum pressure. 


P Strain Indicator. This reads the strain from the pressure 
transducer. It is calibrated for the pressure range of the transducer, 
making it possible to read confining pressures. The strain indicator 


used here is a SR4 Strain Indicator No. 1721A. 
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UP Electricals liming init. Thismre consists of an electrical 
network of resistors, capacitors and inductance coils, which act as a 
pulsating unit. By adjusting the properties of the components of the 
network, the frequency and duration of loading can be controlled as 


desired. Minimum on-load time : 0.25 seconds. 


R Ultra-Violet Recorder. This is a six channel unit. The 
movement of the ultra-violet beam is controlled by the movement of the 
galvanometer. There are six beams each corresponding to one LVDT or 
Ones iGagscell on onewGL che triaxial ‘cells. ~The’ recording paper is a 
treated paper which is sensitive to ultra-violet light. The speed of 
the recording paper can be controlled. The one used here is manufac- 


tured by SE Laboratories, England, model : UV Recorder 3006. 


S Galvanometer and Signal Conditioner Unit. This is also a six 
channel unit which takes in signals from the load cells and LVDTs, con- 
ditions them for the recorder and then passes them on to the ultra- 
violet recorder. The unit used here is : B&F Instruments, Inc. - 


Transducer Conditioner - model 6-101B1-10-200. 


T Galvanometer Switches. These switches switch on the galvano- 
meter to the various channels. By switching off a particular channel, 
the galvanometer in that channel and consequently the operation of the 
ultra-violet beam for that channel, is cut off. The panel used here 


is : B&F Instruments, Inc. - Meter and Switch Panel model SA364. 


U Channel Selector Switch. This switches on the desired 
channel of the six channels so as to read the excitation voltage for 


the load cell or the LVDT. While it is switched on to a particular 
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channel, the load cell or the LVDT in that channel can be balanced to 


their null point. 


Solenoid Valve (Not shown in the figure). This valve controls 
the flow of pressurized air from the air-pressure reservoir to the air- 
pressure cylinder and piston unit. It also acts as an exhaust for the 
air from the air-pressure cylinder and piston unit. This is a solenoid- 
operated three-way valve, the opening and closing of which is done 
electrically and is controlled by an electrical timing unit. The ones 


used here vare : Hannifin =~model;)CCJ<1-37. Range : 0 to 120-psi. 


Air-Pressure Gauges (Not shown in the figure). These gauges 
are installed after the air-pressure regulators and indicate the air 
pressure being supplied by that particular pipe. These are inaccurate 


at low pressures (0 - 10 psi). Ranges: 0 to 60 psi and 0 to 160 psi. 


FIGURE A.2 shows the principal features of the triaxial cell 
along with a typical specimen set up within it. FIGURE A.2 should be 


referred to for corresponding numbers. 


1 stainless-steel loading ram, 0.75 inch in diameter 

2 air release valve 

3 rubber O-rings 

4 0.5 inch diameter stainless-steel tie bars, 6 equally spaced 


5 end cover plates (4 inch diameter x 0.25 inch thick metallic 


discs, for 8 inch high samples) 
6 drainage and pore pressure connection 
7. base of the triaxial cell 


8 connection to pressure supply (for confining pressure) 
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FIGURE A.2 - THE TRIAXIAL CELL (FOR 4 INCH DIAMETER SAMPLES) 
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A.2. 


A.2 


te 


10 


11 


Te 


13 


14 


15 


16 


A9 


4 inch diameter sample enclosed in two rubber membranes 


perspex cylinder, 7.75 inches outer diameter, 7.0 inches 


inner diameter and 11.375 inches long 


reinforcing bands of resin bonded glass fiber, 1 inch wide 
and 0.1 inch thick 


loading cap 

collar and pressure seal 
perspex o11 retainer 
top. capwor triaxial cell 


polythene tube connection to oil supply. 


Operational Procedure 


Setting up the Sample in the Triaxial Cell 


1) 


2) 


3) 


4) 


5) 


6) 


7) 


Place.a.4; inch diameter, cover plate.on; the.seat.in,the 


base of the triaxial cell. 


Place the sample (4 inch diameter x 8 inch high) on the 


cover plate. 
Place a 4 inch diameter cover plate on top of the sample. 


Place the loading cap, with the slot facing up, on the 


cover plate. 


Enclose the sample with two 4 inch diameter rubber 


membranes. 
Align the assembly for verticality. 


Check and see that there are no kinks at the junctions of, 
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8) 


9) 


10) 


14) 


12) 


13) 


14) 


Al10 
a) the triaxial base seat and the cover plate, 
b) the cover plates and the sample, 


c) the cover plate and the loading cap. 


Fix two O-rings, over the rubber membranes to, 
a) the triaxial base seat at the bottom, 


b) the loading cap at the top. 


Fold the excess of rubber membranes around the O-rings, 


neatly over the O-rings. 


Mount the perspex encasing on the base of the triaxial cell, 


and give the mounting screws a few turns. 


Check the eccentricity of the loading cap. This is done by 
lowering the loading ram into the slot of the loading cap 
and rotating the ram. If there is no lateral movement of 
the loading cap then there is no eccentricity. But if 
there is a lateral movement of the loading cap, then 
eccentricity is present. In-this case the perspex encasing 
is removed and the sample and the loading cap are realigned. 
iUsmpLOcessat Sevepeatcd Cillethere ds, NO eccentraci cy. 


present. 
The mounting screws are tightened. 


Close all the connection valves in the base of the triaxial 


cell 


Connect the distilled water supply pipe to the lateral 


pressure supply connection in the triaxial cell. 
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15) 


16) 


17) 


18) 


19) 


20) 


21) 


22) 


All 


Open the lateral pressure supply connection valve and the 


air release valve and supply distilled water to the cell. 


Allow the distilled water to rise to the level of the 


perspex oil retainer. 
Shut the distilled water supply. 


Pump 0il through the polythene tube to the perspex oil 


retainer until the oil surrounds the loading ram. 
Resume the distilled water supply and fill the cell. 


Shut the lateral pressure supply connection valve and the 


air release valve. 
Disconnect the distilled water supply pipe. 


The triaxial cell is now ready to be placed in the loading 


bay. 


Preparation of Bays for Testing 


1) 


2) 


3) 


Close the bleeding valve of the air-pressure reservoir. 


Fill the air-pressure reservoir tank with compressed air 
from the air-pressure line of the building, by opening 


the air-pressure valve. 


Switch on the, 


a) Galvanometer and Signal Conditioner unit. 


(Transducer Conditioner unit.) 
b)ebiectricalo timing unic. 


c) Ultra-Violet recorder. 
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4) 


5) 


6) 


7) 


Al2 


Push the 'Strike' button of the ultra-violet recorder after 


it is switched on. 
Set the paper speed of the ultra-violet recorder. 


Determine the test conditions for each bay (i.e., deviator 
stress, confining pressure, the frequency and duration of 


the load). 


Set the electrical timing unit at the decided on-load and 
off-load times. (The sum of on-load and off-load times 
gives the duration of one cycle, and the on-load time 

gives the duration of the load.) In this study, the on-load 
time was 0.25 seconds and the duration of one cycle was 3 
seconds. FIGURE A.3 shows a typical load-deformation trace 
obtained from the ultra-violet recorder for the above 


conditions. 


FIGURE A.3 - REPRESENTATION OF STRIP CHART FROM ULTRA-VIOLET RECORDER 
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a) Setting up Bay No. 1 


1) 


2) 


3) 


4) 


5) 


6) 


7) 


8) 


9) 


10) 


11) 


12) 


Place a block of wood with a loading cap (dummy sample) in 


the position for the triaxial cell in bay no. l. 


Adjust the loading yoke so that the load cell makes 


proper contact with the loading cap. 


The channel selector switch is switched to channel no. 2. 


(Channel nos. 4 and 6 for bay nos. 2 and 3, respectively.) 
Switch on the galvanometer switch no. 2. 
Check the set sensitivity (6 volts) for disturbances. 


Balance the galvanometer to zero, by adjusting the balance 


knob=nos” 2: 
Switch off the channel selector switch. 


Bring the marker ray for the load cell to the base line by 
turning the galvanometer head no. 2 (located in the ultra- 
violet recorder) with a non-metallic screw driver supplied 


with the unit. 
Switch on the bay operation switch for bay no. l. 


Gradually open the valve (by turning clockwise), which 


controls the air pressure for the deviator stress. 


The magnitude of the load is read from the movement of the 


marker ray. (Set calibration 2:1 Div. = 25 lbs.) 


After the desired load is reached, lock the air-pressure 


valve with the lock nut. 
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13) 


14) 


15) 


16) 


17) 


18) 


19) 


20) 


ZL) 


22) 


23) 


24) 


25) 


26) 


Al4 
Switch off the bay operation switch for bay no. l. 
Leave the galvanometer switch no. 2 on. 
Remove the dummy sample from bay no. l. 
Place the prepared triaxial cell in bay no. |. 


Adjust the loading yoke so that the load cell sits properly 


on the loading ram. 


Lightly clamp the LVDT to the clamp bar, which is already 


clamped to the loading ram of the triaxial cell. 


Screw in the moving rod of the LVDT (provided with the 
screw), into the: threaded Slot. in the top cap of. the 


Criaxtaiccel LL. 


Switch on the galvanometer switch no. 1. (Channel nos. 3 


anges fOr DayenOs.. 2 and 5, TESpectively.) 
Switch on the channel selector switch to channel no. 1. 
Check the set sensitivity (6 volts) for disturbances. 


Balance the galvanometer to zero, by loosening the LVDT 


clamp and moving the LVDT up or down. 
Clamp the LVDT firmly after the galvanometer is balanced. 
Switch off the channel selector switch. 


Bring the marker ray for the LVDT to the base line by 
turning the galvanometer head no. 1 (located in the ultra- 
violet recorder) with a non-metallic screw driver supplied 


with the unit. 
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2) 


28) 


29) 


30) 


31) 


32) 


33) 


34) 


0%) 


36) 


A15 
Leave the galvanometer switch no. 1 on. 
Set the three-way transducer switch to bay no. l. 


Connect the lateral pressure supply line and the pressure 
recording line for the transducer, to the lateral pressure 
supply connection and to the pore-pressure measurement 


connection, respectively, in the base of the triaxial cell. 


Open the valve in the connection, connecting the cell 


pressure to the transducer. 
Read the strain indicator for zero cell pressure. 


Convert the desired cell pressure to be used, to the strain 
indicator’ reading. (Set sensitivity : 1 psi = 33.2 div. of 


strain indicator for transducer no. 28749.) 


Add the converted strain indicator reading to the reading 
FOtrezeroecctiapressure sand set, the strain indicator dial 


to that value. 


Slightly open the valve controlling the air pressure for 
confining pressure (by turning clockwise), and then open 
the lateral pressure supply connection valve in the base of 
the triaxial cell. (This is done to avoid water from the 


triaxial cell coming into the air-pressure line.) 


Continue increasing the air pressure supply till the meter 


in the strain indicator is balanced to zero. 


Lock the air-pressure valve with the lock nut. 
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37) 


38) 


39) 


40) 


41) 


42) 


Al6 
Set the counter to zero. 


Adjust the guides on the loading yoke (by turning them on 


the thread) so that they come within the guide frame. 


Determine the anticipated maximum deflection the sample is 


likely to undergo. 


Adjust the gap between the safety microswitch and the steel 


bar according to the determination in the previous step. 
Bay no. 1 is ready for operation. 


To start operation, switch the paper movement switch in the 
ultra-violet recorder to 'Drive' and then switch on the bay 


operation switch for bay no. l. 


b) Setting up Bay Nos. 2 and 3 


The procedure for bay no. 1 is repeated for bay no. 2 and bay 


Mle, SN 


A.3. Precautions 


1) 


2) 


The sensitivity control knobs of the galvanometer and 


Signal conditioner unit should not be touched at all. 


Once the test is started, the LVDT should not be handled 
unless the marker ray is about to go off the recording 
paper“and requires resetting. In that case, the reading 
of the LVDT marker ray should be taken and then the LVDT 
should be reset so that the marker ray comes to the de- 


Sired base line. It must be noted that on this new base 
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3) 


4) 


5) 


6) 


7) 


8) 


Al7 
line, the value of the last LVDT reading must be added to 
the reset base line (which for all practical purposes is 


taken as zero, for the new operation). 


Adjust the galvanometer heads with a non-metallic (non- 
Magnetic) screw driver supplied with the unit, and not with 


any metallic screw driver. 


Never remove the recorder cover without having the recorder 
turned off because, 


a) The fan is very close to the cover and may be damaged 


if the cover is removed while it is running. 

b) Light damage to the eye would be severe, if looked 
directly at the bulb. 

Dust cover should be placed over the ultra-violet recorder 


when it is not in use. 


Always switch off the various power switches in individual 


units, before disconnecting the main power plug. 


Precaution should be taken never to apply a negative 


pressure to the pressure transducer. 


All electronic equipments and parts functioning electroni- 


cally must be handled delicately. 


A.4. Electrical Connections and Calibration Curves 


The remaining pages of this appendix give the schematic diagrams 


of the electrical connections for the measuring and recording equipment. 


The calibration curves for the load cells and the pressure transducer 


are also included. 
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The uncertainty in the laboratory experimental results were 
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calculated taking into account the uncertainty in each measuring device, 


by the root sum square method of combination*. The equation for calcu- 


lating the resilient modulus is 


Oo 
M a A/1 eoererereeeo reese eee ere sees eee ee eee eeeees Ji\es it 
in which 
M = resilient modulus, 
o = deviator stress, 
A = recoverable deformation, and 
1 = length of sample. 
The uncertainty in resilient modulus is given by 
1/2 
W. cs OM , oe OMe ee CP er Z Nee? 
M oO O ol 1 0A A : 
in which 
ie uncertainty in resilient modulus, 
We uncertainty in the measurement of deviator stress, 
ae uncertainty in the measurement of length of sample, and 
Wer uncertainty in the measurement of recoverable deformation. 


* J. P. Holman, “Experimental Methods for Engineers", 2nd ed. 
McGraw-Hill, 1971. 
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Sample Calculation 


Problem@paramecers are; §o = 8 psi,-l)= 8 ins., A’ = 0.0075 ins., 
Were 7 Oe DSis CAP POO; OF eins. 
and Waae * 0.5% = 0.00004 ins. 


From Equation A.1 and problem parameters 


3M poll 

ae 1066 

0M ej 

ORK eat ee 1066 

OMyy ig 6 Guknlyis 1138000 
oA Re 


Substituting in Equation A.2 


= 
i} 


1/2 
a0se x 0.2)2 + (1066 x 0.01)2 + (1138000 x 0.00008) } 


45242 si 


From Equation A.1 and problem parameters 


M = 8550 psi 
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APPENDIX B 


AUTOMATIC KNEADING COMPACTOR 


B.1. Description and Features 
An Electronic - Hydraulic Kneading Compactor (model CS 1000 


manufactured by Cox and Sons of Sacramento, California) has been set up 
in the laboratory at The University of Alberta. Although the compactor 
is basically designed for requirements specified in California Division 
of Highways Test Methods 301 and 304, the specimen size normally being 
4 inch diameter by 2.5 inch high, it can be used for other compaction 


procedures and specimen heights up to a maximum of 8 inches. 


The apparatus is capable of imparting nearly identical proper- 
ties to a number of samples which are compacted separately in the 
machine under the same initial conditions (such as moisture content, 


6teu) 


The tamper foot may be raised or lowered through a maximum 
height of 9 inches with individual adjustments of : down stroke rate, 


up stroke rate, up stroke return distance and dwell in down position. 


The rotating table is hydraulically actuated and electronically 
timed to the tamper foot. The table rotation can be varied from 6° to 
72° increments for each 360° of table rotation, by adjusting the ten- 


turn control dial for table rotation. 


The automatic feeder assembly is an endless belt type, activated 


by a hydraulic actuator and electronically timed to the tamper foot. 
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B2 
The loader-conveyor belt feed may be stopped at any time during a load 
cycle without stopping compaction of the already deposited material. 
The count will continue but no additional material will be fed from the 


belt during the interruption. 


The return stroke, ram dwell, table rotation and loader stroke 
can be changed with time to meet the requirements of a compaction 
procedure. Timing periods are set with precision ten-turn controls and 


timing of each mode is independent of all other adjustments. 


The compaction control mode is divided into three functions, 


a) Load A, b) Compact, and c) Load B. 


The automatic load cycle A permits the specimen to be loaded 
into the mould in up to 30 equal partitions, with compaction tamps 
interposed between each incremental addition to the mould. At the con- 
clusion of the loading process, the machine automatically terminates 


operation. 


The automatic load cycle B uniformly loads the specimen into the 
mould in two volumes of 10 partitions each. After the first half is 
evenly spread around the mould, the automatic sequence provides 10 
compacting tamps of the ram. The second volume is similarly loaded, 
followed by 10 additional compacting tamps, after which the machine 


automatically terminates operation. 


The compact mode will continuously tamp the specimen until 


terminated automatically by the predetermining counter or operator. 
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The Load A and Load B modes are primarily meant for the standard 
California test methods, but could be used for other purposes as well, 
if found suitable. Normally the Compact mode is used for different 


compaction procedures. 


FIGURE B.1 shows the essential features of the automatic 
kneading compactor and it should be referred to for the corresponding 


alphabetic letters while reading the functions of the various components. 


Bee. 1. Functrons 
A Ram Pressure Gauge. Measures the pressure on the specimen 
directly in pounds per square inch. Total force equals 3.2 foot area 


multiplied times the gauge reading. 


B Predetermining Counter. This counter is used to preset the 
number of compactive tamps to be applied to the specimen. The counter 
must be set at zero to begin fabricating a specimen. To preset, 


1) Release the locking mechanism of the lid by pressing 


pushbutton. 
2) Open the lid and set figure wheels. 
3) Close the lid. 
4) The desired number is now visible in the window. 
5) The figures must be in alignment. 


6) Resetting to zero for repetition or change of pre-selected 


number is done by pressing the pushbutton. 


7) CAUTION: DO NOT ATTEMPT TO PRESET OR RESET COUNTER 
WITH MACHINE RUNNING. 


8) Ram will not operate until the counter is reset. 
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FIGURE B.1 - AUTOMATIC KNEADING COMPACTOR 
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C Counter Indicating Lamp - Green. An intermittent light that 
flashes each time the ram advances toward the specimen. Indicates that 


switch D is in On position. 


D Counter Off-On Switch. Must be in On position for the 
counter to operate automatically. Off position permits ram operation 


without counter control. 


E Return Stroke. This ten-turn dial controls the time interval 


between strokes of the Ram. 


F Ram Dwell. This ten-turn dial controls the amount of time 


the Ram applies pressure to the specimen. 


G Master Indicating Lamp - Red. A light that indicates 


whenever the master switch is in the On position. 


H Master Switch. Controls the electrical power to the machine. 
Prior to turning off the master switch, ALWAYS place switch L in 


Ram-Up position. 


I Loader Pulse Button. Used to advance the conveyor belt of 
the loader manually. Load pulses may be interposed by the operator at 
will by means of this manual button. Loader switch J must be in the On 


position. 


J Loader Off-On Switch. Must be in On position for the 


control system to activate the conveyor belt. 


K Reset Button. Initiates the start of all control modes. 
Example: Place switch L in Compact, switch M in any desired mode, and 


push reset button. 
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L Ram-Up - Compact Switch. Ram-up position retracts the ram. 


For SAFETY always use this position when placing hands under the Ram, 


and before turning off the master switch. 


M Control Mode. This control is divided into three functions 
Load A, Compact, and Load B. To initiate any of the above modes, place 


switch L in Compact position and push reset button K. 


N Ram Pressure. This knob is used to control the amount of 
compactive pressure applied to the specimen in pounds per square inch. 
Turning this knob clockwise will increase the pressure and turning the 


knob counterclockwise will decrease the pressure. 


O Heater. This is the tamper foot heat control for Bituminous 


Mix specimens. Turn clockwise to increase temperature. 


P Heater Indicating Lamp - Amber. This light indicates when 


the heater is on. 
Q Heater Off-On Switch. 


R Table Rotation. This ten-turn dial controls the amount of 


movement of the table between tamps of the Ram. 


S Loader Stroke. This ten-turn dial is used to control the 


length of movement of the conveyor in the loader trough. 


T Down Rate (Ram velocity). This knob is used to control the 
downward velocity of the Ram. Turn clockwise to decrease, counter- 


clockwise to increase rate. 
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U Up Rate (Ram velocity). This knob is used to control the up- 
ward velocity of the Ram. Turn clockwise to decrease, counterclockwise 


to increase rate. 


V Automatic Conveyor. This device is preset by electronic 
controls to load in the proper sequence. 

Range of belt movement : 0.12 to 3.00 inches per pulse. 

Capacity of Conveyor Belt : 1,500 grams. 

W Height Indicator Device. This device is not used to 
actually measure the height of a specimen, but rather to give the 
operator an approximate idea of the specimen height so that he can be 


as close as possible to the proper height with the first specimen. 


X Ram. The ram transmits the hydraulic pressure of the 


system to the specimen. 


Y Table. A rotating platen upon which the specimen is 
compacted. This table rotates when the ram is retracting after having 


applied its compactive force. 


To change setting of the ten-turn dial controls EH, F, R and S, 
1) Release dial lock by pushing counterclockwise. 
2) Choose desired number setting. 


3) Reset dial lock. 
B.2. Operational Procedure 
1) Decide on the compaction procedure. 


2) Switch on the Master switch (H), and allow the machine to 


warm up for about 5 minutes. 
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3) Decide on the settings for Return stroke, Ram dwell, Table 
rotation, and Loader stroke and set the respective ten-turn 


dials accordingly. 


4) Up rate and Down rate settings for the standard California 
test methods are adequate for normal use. If a change is 


desired, they should be set accordingly. 
5) Set the mould on the Table (Y). 


6) If the mould is to be loaded by conveyor belt feed, place 
material on the conveyor belt and spread it to a uniform 
level. Swing the conveyor assembly to feeding position and 


switch on the loader switch (J). 


7) If the mould is to be loaded manually, place the material 
in the mould and spread it out evenly. Step 6 is omitted 


fox thisecase. 


8) Select the Control mode (M) to one of the three functions, 
as desired. Load A and Load B cycles function automatically 


for feeding and compacting, as explained earlier. 


9) If the mode selected is Compact and the number of tamps is 
to be controlled by the predetermining counter, set the 
desired number of tamps in the predetermining counter (B) 


and switch on the counter switch. 
10) Place Ram-up-Compact switch (L) in Compact position. 


11) Press the Reset button (K), which starts the compaction and 


all control modes. 
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12) The pressure applied by the tamper foot is read in Ram 
Pressure Gauge (A). The Ram Pressure knob (N) is manipu- 


lated till the desired foot pressure is reached. 


13) After completion of the cycle the Ram retracts and machine 
terminates operation, the Ram-up-Compact switch (L) should 


be placed in the Ram-up position. 


14) If another layer of material is to be added manually and 
compacted, steps 9, 10, 11, and 13 are repeated after the 


material is placed in the mould. 


15) After the sample has been fabricated and the machine is no 
longer required for operation, place Ram-up-Compact switch 
(L) in Ram-up position and switch off all the switches, 


the Master switch (H) being switched off last. 


B.3. Procedure for Fabricating 4 inch Diameter x 8 inch High Specimens 
The procedure for fabricating 4 inch diameter x 8 inch high 
samples was different from that of smaller samples. The normal working 
of the compactor was restricted due to the large sample height 
(8 inches). The loader-conveyor belt feed could not be used. A 
clearance of the cross-arm of the Ram above the top of the mould was 
only obtained after a 1.5 inch layer of compacted material was inside 
the mould. It was evident from this difficulty that the sample would 
have to be compacted in five layers or less, to have the machine 


working properly. 
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After several trials it was decided that the material required 
for preparing one sample be weighed into 5 equal parts. Each part was 


fed manually into the mould and given 25 tamps of the tamper foot. 


During the trials, it was observed that the foot pressure 
greatly influenced the foot penetration when the moisture content of 
the soil was wet of optimum moisture content. Hence different foot 
pressures were used for calibration, depending on the moisture content 


of the soil, and the trial results are summarized in TABLE B-l. 


A series of compaction tests were also run with a different 
soil at various water contents, for which the ASTM D698-70 Test Curve 
was available. The mould used was 4 inches in diameter and 4.5 inches 
high, similar to the one specified for the Standard ASTM D698-70 Test. 
Samples were compacted in three layers with the machine settings 
indicated in column (2) of TABLE B-l. The results obtained were com- 
pared with those obtained from the Standard ASTM D698-70 Test, and are 
shown in FIGURE B.2. 
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FIGURE B.2 - COMPARISON BETWEEN ASTM D698-70 AND KNEADING COMPACTION TESTS 


% 


bot iupe ‘yrsism oft tert behiosh esw Ff Shaelrd [aveved TOTTA 
ise4 .@teeq [aupe 2 otk Bedgtow oc oiqnee eno anineqorq rot 


2S mavig bos Sivom oft otni Ltaunsm bok 


a, 7 


n =" sola " f pune t ofr . e - ar fue 4 » ofl ~ i = 
3rq foot oft studs bevrerda 28w ti, ebeirs ods gaimd 7 
‘ — . ‘ a ~ 
y ik Hoke MIDE 20Q1 art b someul> nt Uissexy 
o? F (bh oogoH .tnetnod emrielom mumitgqo te 3ew aeaw [foe * 
=e 
ion oft 7 gnibssyesh .moisetdilas tot boaw erew aeTuezetg 


% 
| oe 
ni bovivamme sta 2p iweer Leitt edd bae .ltoz off ie 
ca 7 


yor 
aie Mee ‘ + ¥ 4 es ' ~ we 3 t r; 
4 | ] 7 eae. | | -Je03 AWTS. OMmoOos a asitec A : 7 
; i 
~*~ ry 
i : 
PT. T rp set dinibelas at Straten” « bee bey Vs o 
{ \ ian Shs tees f! IG A © eINSInesd ) im BO 2uoirs av J& heVe 
t ra ' 
. 29xk9 bh apw beso bluoam oaT- .oldaf. 
win 
‘7 r — 7 ~ ‘ 
; $2 ent 03 aq 10 of? ot tel imte 


4 = 


moy stew beristdo etives: edf .f-@ GIZAT Zo (5) maylos at 


©. qaAUOr4 ai m 


~~ 


~~ 
re 
; 
g 
re 
= 


~ 
ae 
‘ 
= 
eo 
a2] 
, - 
meTeyr yp 
is se 
i 
¢ iu 


N\tr. 


, 
c x 
" 

: ai 


° 8800 MTGA 
werd bb 20T 


Py. at yf “sf 


; cn gnibeoe 


Bll 


qua quo) aanjsyoR wnsy3dg - ONO 
BBTMYIOTIAGIUNOD - MDD : ALON 


*uotTqoedwo. 
JO aT9h>d AayQOUP OF 
peaqoafqns sem o{ dues 
ay. s8uy3jIas AaYyIO 3y3 
Buyi1azTe InoyI}TmM pues 
tsd Q¢ 02 paonped sen 
@anssaad 300J ey J pe3I9 
-edwod sem zohel TeUuys 
ay2 aeqagy ‘*seyout 
6Z7°0 sem uot Ierjaued 
Qoog *saakeT ¢ ut 

pe qoeduo> sem aT dues 


*mOT Sea uoTIeIJaUed 
yooy *saakel ¢ Ut 
paaedwod sem e,dwes 


*Juaquod aianjstow 
% O'St 38 °30d g°60T 


*qyuaquos sianqysjou 
Ch ata GRD Daiexel Pa2reigt 


01 - 869 d WISV PrBpueds 
aya Aq pauteygo Aqysuap Aig 


*quaqjuod sianystouw 
% Z°ST ae ‘3d LOTT 


* quo quo5 3ainjsjyou 
%1°ST ae *3>d €°ZTT 


*s8uy 3396 aAcgeB 
au? YITM 3IOJDedWOD ZuTpEaUuy 
aya Aq pauteqgo Aqysuaep 41g 


24907 
mod suana ¢°¢ 
0z0 
Oct 


peyxd07T 
MOO sUInd ¢’°E 
ost 
Ov 
OOT 
eTPI 
4 
38d O00T 


P2490 
MOD suing ¢°E 

OST 
Ov 
001 
eTPI 
SZ 
ysd OST 


aqey umog 


Treg wey 
ayOIIS UINIDY 


uot zeIYy sTqQeL 
(s) 3y019S$ JepeoT 
Jaqunoy) Zuyuywieyeperg 


ainsseig 2004 


(1) 
9WO jo Aap io 
OWO 3B 6TFOS 


jaanjoejnuew ay3 Aq pattddne se 


UOTIBIGILeD eTUIOJTTeD paepurys 


OHO JO 394 BT TOS 


YOLIVdWOD ONIGWANY JHL ONISN “SNSWII4dS HOIH HONI 8 A@ YALAWVIG HONI V7 
ONILVOTUdV4 YO4 SISAL NOILIVdWOD WWIul JO AYWWWNS 


L-d@ 31avl 


; agisewiled elemilis) besberast t 7H) fo gov ad tot | We ge aliat. 
; tetutsetimes eds Ys ball gque cas | , ; MS to ¢xk vo : 
. i : : a 
; . : ? 5 
| ‘ (e ft) ; (t} : ; 
; < i - << 4 a - 
jo ne ee ee of in eee Se so ee a aprenerion oe -hos 
: ; ; ; 
* 4 2eq «400% - /$eq G2i () 
: $ ‘ : 
. a wo | 
; « : ill nbs : aX 
i - : P (23 
| - : ; 
: os Gai oe | 48> 
: “eG f cei 
; 
| 809. eocwz 2.€ We eowe 2.6 WO9 unzes 2.£ ; 
; Bais ; : iw 
bs: ; j : F ; 

ee Sa poo a ath Set es << me seer obama: 
: ; : 
i Le (i baqiecdo ytiassh = 

> ; / ; # ud f.ci-3 °.81i | s47 Shiv ms2eqias soibese 

: oa 7 i923 saget H -tpoisje8 Bee 

; ; 
- : on x 2 6 ” =f * 

- - is Gq <«346 > o 
: ; 
; ' » 4 
; P ; 
F é @ 
: : 
‘ : 
: f . 
: t 

: 


o- 


“DRITADIASAT AOW- great norTOAa MO JAIA 
AOTIAIMOD “SVIGAINM 


Se 


a 


JHT AKICU evaMi9392 Ho TH HOMI 8 


oe 7 


T 0 val weer 


B12 


B.4. Recommended Test Procedure for Obtaining Standard ASTM D698-70 
Densities in 4 inch Diameter by 8 inch High Specimens 


The specimen is fabricated in 5 layers with 25 tamps per layer 


and machine settings set at, 


Predetermining Counter 
Loader Stroke 
Table Rotation 
Return Stroke 
Ram Dwell 
Down Rate 
Up Rate 


Foot Pressure 


25 
Idle 
100 
140 
150 
53.5 turns CCW 
Locked 


Variable depending on the 


moisture content of the soil. 


(150 psi at OMC and 100 psi 
at about 3% wet of OMC). 
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APPENDIX C 
WANG 720 C COMPUTER PROGRAM FOR PROCESSING 


REPEATED LOADING TRIAXIAL TEST DATA 


C.1. Abstract 
The program computes the total, residual, and resilient strain 
in percentile form and the resilient modulus in psi for a particular 


number of load repetition, and prints the computed results. 


C.2. Program Initiation 

The program is in the first block in Side 2 of Mp Program-Tape. 
In case the Mp Program-Tape is not available, the listing of the pro- 
gram given in TABLE C-1 should be keyed into the Wang 720 C Computer 
and the program recorded in a tape. Program Loading Steps; 

1) Primes 

2) Load Program (Block 1 in Side 2 of Mp Program-Tape). 

3) Verify Program 3159. 


Details are given in section C.4 


C.3. Program Description 


Usually the total strain, residual strain, resilient strain and 
the resilient modulus are plotted against the logarithm of the number of 
load repetitions, and this program computes these values and prints 
them in a format from which it is fairly easy to plot these graphs. The 
number of load repetitions is programmed in a logarithmic scale such 
that the value is automatically incremented to the next number of 


repetition and is displayed in the y-register after the machine 
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completes the computation of the output parameters for a particular 


number of load repetition. 


The input for the program are the sample number, the deviator 
stress, the confining pressure, the LVDT sensitivity, the sample length, 
the total strain and the residual strain. The LVDT sensitivity is in 
the form of one division in the strip chart representing a certain 
length (in inches). As such, the total and residual strains from the 
strip chart are read in the form of number of divisions and are the 
input parameters for the program. It is from all these inputs that the 
program computes the respective results. A representation of the strip 
chart from the ultra-violet recorder is shown in FIGURE A.3 in 


APPENDIX A. 


C.4. Operating Procedure 
1) Load the program tape, depress REWIND, depress TAPE READY, 


depress RUN, PRIME, LOAD PROGRAM. 
2) PRIME, VERIFY PROGRAM; x-register should read 3159. 
5) = Printer ON> AUTO. 


4) PRIME; key special function 00 01 
Printer will respond with : 


eoeoeeeree eee eee eee eee eee eee ee eee 


5) Displaveeaws Lo- ey 
2- X 


Enter data: First part of sample number in y-register, 


Second part of sample number in x-register. 
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Key GO 
Printer will respond with : 
SAMPLE NO.L Y/X 

Display 4 =] 

4-X 
Enter data: og in y-register, 

Oz in x-register. 
Key GO 


Printemowi lJ) respondswith. : 
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Note: -d and 3 must be prefixed with o, and since the 
printer does not have that character it must be 
Wratten later. 

Display oq) bode 

6 - X 
Enter data: LVDT sensitivity in inches in y-register. 


Length of sample in inches in x-register. 


Key GO 
Printer eerlaLe: 
Display ee! Ee ¢ (displays the number of repetition) 


T= X 


Enter data: Total strain reading in x-register 


(number of divisions in strip chart). 
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Key GO 
Printer ;/etdie. 
Display Sk sey (displays the number of repetition) 


8 - X 


Enter data: Residual strain reading in x-register 


(number of divisions in strip chart). 


Key GO 

Printer meeWitdoprint the caletlated$totalestrain, 
residual strain, resilient strain, and Mp 
corresponding to the number of repetition of l. 
See TABLE C-2 for printer outputs. 

Display. dan2ia J. 


Vi eS 


Repeatestenss il, 12,13, andml4e 


if at stepais at as realized that wrong data, was,entered 
for total strain in step 11, press keys SEARCH, 7 and then 
enter the correct data and proceed to step 12 (the program 


back#éteps*totstep unl) . 


The number of repetition is incremented in a logarithmic 
fashion after a set of data have been entered and computed. 
it is incremented in intervals of 1 from 1 to 10, in 
intervals of 10 from 10 to 100, in intervals of 100 from 
100 to 1,000 and so on. If the display in the y-register 
at step 11 shows a repetition for which there is no data, 
press keys SEARCH, SEARCH and the number of repetition will 


be incremented to the next value. 
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If the number of repetition is different from the 
normal programmed logarithmic scale, enter the actual 
number of repetition in y-register before entering the total 


strain reading in the x-register at step 11. 


Example 
Display atestep 11: 700 - Y 
7 -X 
Programmed sequence Actual sequence 
of No. of Rep. of No. of Rep. 
600 600 
700 900 
800 950 
900 1000 
1000 


1) Entex 900 in y-register and do steps 11 to 14. 


Display : 800 - Y (Note that the number of repeti- 
7X tion is still in the programmed 
sequence) 


2) }Enter 950 in y-register and do steps 11 to 14. 


Display =: 90054 
7-X 


3) Key SEARCH, SEARCH 


Display : 1000 - Y 
7 - X 


Proceed normally from step 11. 


livates Cepsto it 1s found) tatethne nexcese tr OLecacdator 
total strain and residual strain is the same as those just 
entered, press keys SEARCH, GO and the printer will print 


the same results again for the number of repetition 
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C6 
currently displayed in the y-register. However, if the 
number of repetition displayed in the y-register at this 
stage is other than what is required, then it should first 


be changed according to note b). 


Example 
Data 
Sample No. : L 6/1 
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LVDT sensitivity : 1 div. = 0.001 inches 
Length of sample : 6.135 inches 


From the strip chart of the ultra-violet recorder : 
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Do the steps described in section C.4 (operating procedure). 
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APPENDIX D 


STIFFNESS OF ASPHALT MIXTURE BY INDIRECT METHOD 


D.1. Stiffness of Mixture - Loading Time Relationship 

The procedure used for the determination of stiffness of 
mixture - loading time relationship for the asphalt concrete used in 
the construction of Highway 15-A-1, is based on the procedure suggested 
by McLeod (1969). A sample calculation for determining the relationship 
at 60°F is included, which constitutes the master curve shown in 
FIGURE 5.1. The same relationship was determined for other temperatures 
and using the shift factor concept as described by Finn (1967), these 
curves were transferred to the master curve. The shift factor - temper- 


ature relationship is shown in FIGURE 5.2 in the text. 


D.2. Determination of Master Curve and Shift Factor 


eco Master, Curve 
The recovered properties of asphalt concrete is presented in 
TABLE III-2 in the text, from which the relevant properties required 


for the determination of the master curve are: 


Penetration at 77°F (25°C) ; 169 
Viscosity at 140°F (60°C) : 572 poises 
Percent air voids OP EE) 
Percent voids in mineral aggregate (VMA) 2416.5 
Volume concentration of aggregate (C,) = 0286 
Corrected volume concentration of aggregate (C') : 0.84 
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D2 
Procedure: Refer to McLeod (1969) page 172. 


Step 1: From recovered properties, 
Penetration at 77°F : 169 


Viseasity at 140 F-92572 poses, 


Step 2: From these values, 


Penetration Index ree etal aa 0 Bs 


Step 3: For penetration at 77°F = 169 
and penetration index = -1.0 
Temperature difference obtained = 14°C 


Hence Base temperature = 25 + 14 = 39°C. 


Step 4: The curve is to be determined for a service temperature 
ofs60Fa(l5.6)C) Sewhich is 23.45C.(39.0 - 15.6 = 23.4) 
below the Base temperature. 

For a loading time of 0.01 seconds, stiffness modulus 


of asphalt cement = 200 kg/cm? = 2.84 x 10 psi. 
Step 5: Volume concentration of aggregate, 


_ 100 - VMA Paleo 5 na 
vane VOOR Maltevolds = slO00s-we3.0> 4 wi 


Corrected volume concentration of aggregate according 


to Van Draat's and Sommer's equation, 


C 
aie Ree ae a EO eer ee 
v ~ Trav ~ 1 (0.057 - 0.030) 
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Step 6: Stiffness of asphalt mixture at 60°F, 5.3% asphalt 
content and 5.7% air voids corresponding to a 
stiffness modulus of asphalt cement at 60°F and 
loading time of 0.01 seconds = 2.84 x 10° psi, 


is 300,000 psi. 


Similarly, stiffness of the same mixture at the same temperature 
(60°F) and different loading times were calculated and are tabulated in 


TABLE D-1. 


TABLE D-1 
STIFFNESS OF MIXTURE FOR VARIOUS 
LOADING TIMES AT 60°F 


Loading Time Stiffness of Mixture 
Secs psi 
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These values were plotted (both in logarithmic scale) to give 
a stiffness of mixture - loading time relationship, and yielded the 


curve presented in FIGURE 5.1. 
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Wed. 2. ) ohitt Factor 

The whole procedure was repeated for different temperatures and 
the amount by which the curves were shifted horizontally (termed as 
shift factor) to overlap on the master curve were found and are 


tabulated in TABLE D-2. 


TABLE D-2 
SHIFT FACTORS FOR VARIOUS TEMPERATURES 


Temperature Shite. factor 
°F 


so4axe 105+ 


.00 x 10° 
ppOex LOL 
55 x 101 


These values were plotted (shift factor in logarithmic scale 
and temperature in linear scale) to give a temperature - shift factor 


helationship and yieldedsthe curve presented in FIGURE 5.2. 


Dec. Use Of sche Curves 
An example is illustrated on how these curves are used to 


determine the stiffness of the mixture. 


Problem : Determine the stiffness of the mixture at a 


temperature of 42°F and loading time of 0.1 seconds. 


Solution: Determine the shift factor from FIGURE 5.2. 
For a temperature of 42°F, Shift factor = 8.0. 


Calculate the modified loading time (the loading 
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DS 
time at 60°F for the actual loading time at 42°F) 
from the equation : 


Actual loading time 


modified loading time = Shift Factor 


For actual loading time of 0.1 seconds at 42°F, the 


modified loading time = a 1.25 x 10°72 seconds. 
Enter this loading time in FIGURE 5.1 and find the 
corresponding stiffness from the curve, which is : 


270,000 psi. 


The stiffness of the mixture at 42°F and for a loading time of 


051 seconds 15° 270,000 “ps1. 
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APPENDIX E 


CHEVRON COMPUTER PROGRAM 


Seeeee MAIN ROUTINE - N-LAYER ELASTIC SYSTEM *#4¥*%4# 
DIMENSION RR(29),2Z(20),E(5) /V(5) gHH (4) »H (4) ,AZ (190) ,a(190,5) , 
1 B(190,5), €(190,5),D(190,5) ,AJ (190) ,RJ1(190) ,RIO (190) 
COMMON RR,Z2,E,V,HH,H,hZeA,ByC,yD,AJ,RI1,R250 
COMMON R,Z,AR,NS,N,L, ITN, 2, RSZ,KST,RSR,RTR, ROM, RMU, SF 
COMMON CSZ,CST,CSR,CTR,COM,CMU 
DIMENSION TITLE(18), BZ(100), X(5—¢4,4), SC(4), PM(4,4,4), FM (2,2), 
VITEST (11) 
COMMON TITLE, PSI, NLINE, NOUTP, NTEST, TEST, 
1 ITN4, EGohIE eTZL2 cPRe MPA, ERPeTAPe \TIM; 21, 612, 23, S48: 
2° Tesere atze serie, WA, BJ1,) 890, BZ, ZF, SZ1,°SZ2, PM, SGi, SG2, 
3 PH, PH2, VK2, VKP2, VK4, VKP4, VKK8, X, SC, FM 


INITIALIZE VARIABLES 
NOUTP=0 
ITN = 46 
ITNY = ITN*Y 
READ (5,310,END=62) (TITLE (I) ,I=1, 18) 
FORMAT (18A4) 
IOK=0 
ICALL=0 
READ (5,311) WGT,PSI 


FORMAT (2F12.0) 

READ (5,312) NS, (E(I),V(1), I=1,NS) 
FORMAT (12,5 (F8.0,F6.0)) 

NSSENS shi 

READ (5,313) (HH (I), I=1,N) 


FORMAT (4F6.0) 

IF (NPUN) 7,7,8 

PUNCH 355 

READ (5,301) IR, (RR(I) ,I=1,IR) 
READ (5,301) TZe(22 (yg t= 1,912) 


FORMAT (16, (11F6.0)) 

Ak = SORT (0.31831 ¥*WGT/PSI) 

NLINE = 174NS 

NPAGE = 1 

WRITE(6, 350) (TITLE(I), I = 1,18 ), NPAGE 

FORMAT (1H1//1HO, 18H #*# RRR ERERERERKEKK 61K, 18A4, 1X, SHER EKK KKEEEKEE 
2e***x* 10X, 7H PAGE I3 ) 

WRITE (6,351) WGT, PSI, AR, (I,E(I)/¥(1I) ,HH(1) ,I=1,)) 


FORMAT (1HO, 4OX, 26HTHE PROBLEM PARAMETERS ARE/ 
2 THO; cO Me eI ZHLOTAL LOADS SC, 8X, FI0.2, SH LBS/ 
3 UO PIC 5 USGI Rse IEE SSW RUG ag Bish IS Gr tet eispe7 

4 THO paLO Xe SaLOADSRADIUS. 2707 X— fle 2, soHe INS / sth 

5 Hee 20K, POH LAY CR eel 37 el OH HAS MODULUS gtelio «1 Ol¢ 

6 18H POLS SONS MRALE Olt, me Or Sig mtie AND THICKNESS , FS.1, 

7 GH IN.)) 

WRITE (6,354) NS, E(NS), V(NS) 

FORMAT(1H ,20X,5HLAYER,I3,14H HAS MODULUS ,F8.0, 

1 18H POISSONS) RATIO | Fo.3, 24H AND IS SEMI-INFINITE. ) 
WRITE (6,352) 

FORMAT (1H0,45X,11HS T R ES S,48X,11HS TRATI N/ 

1 1H F 20X HOH FFE RHE HEE EARS EK EKER EERE URE REE REE EEE EEE KEK EK EEE 
QEEKEHEKEKEKE aiGke JOH ERR E RHEE KKK ERE EKER RR KARE HEH KK fi 
Seon R,6X,1HZ,9X,8HVERTICAL, 3X ,1OHTANGENTIAL,7X,6HRADIAL, 8X, 


4 SHSHEAR,OX,4YHBULK,&X,8HDEFLECT ,7X,6HRADIAL,IX,4HBULK/IH ) 
TCON=0. 
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ADJUST LAYER DEPTHS ** 
H (1) =HH (1) 
DO 25 I=2,N 
H (I) =H (I-1) +HH (I) 
IRT=0 
START ON A NEW R *® 
IRI=IRT+1 
IF (IKT-1R) 105,105,10 
R=RR (IRT) 
DONG 1 2Ei=1;12 
DO 31 J=1,N 
TZ = ABS (H(J) - 22(I)) 
LF (T2000 1) 32 532,31 
Z2(I) = -H(J) 
CONTINUE 
WRITE (6,355) 
NLINE = NLINE+1 
FORMAT(1H ) 
CALCULATE THE PARTITION ** 
CALL PART 
CALCULATE THE COEFFICIENTS ** 
DO 125 I=1,1TN4 
P=AZ (I) 
CALL COEE (1) 
IF (R) 115,115,110 
PR = P*R 
CALL BESSEL (0,PR,Y) 
RJIO(I) = Y 
CALL BESSEL (1,P8,Y) 
RJ1(I) = Y 
PA=P*AR 
CALL BESSEL (1,PA,Y) 
AJ (I) =¥ 
CONTINUE 
IZT=0 
START ON A NEW Z ¥* 
IZI=IZT+1 
IF(IZT-IZ) 205,205,100 
Z=ABS (ZZ(12T)) 
IF ( NLINE - 54) 207,206,206 
NPAGE = NPAGE + 1 
NLINE = 8 
WRITE (6,350) (TITLE(I), I = 1,18 ), NPAGE 
WRITE (6,352) 
CONTINUE 
PIND THE LAYER CONTAINING Z ** 
TZZ = 0.0 
DO 210 J1=1,N 
J=NS-J1 
IF (Z~H(J)) 210,215,215 
CONTINUE 
ee | 
GO TO 34 
L=J+1 
IPF (Z2(IZT)) 33,34, 34 
eg 
TZZ°= 1.0 
CONTINUE 
CALL CALCIN (IOK, ICALL,IRT, TCOM) 
IF (122) 36,36,35 
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ZZ (IZT) = -%2 (127) 

I2T = I2T-1 

CONTINUE 

IF (NPUN) 200, 200, 220 

NZP = IzT-1 

NRP = IRT=1 

PUNCH 356, CSZ, CST, CSR, CTR, COM, CMU, N2P, NRP, TITLE(1) 
FORMAT (1X,6E10.3,212, 9X,A6) 

GO 10 200 

STOP 

END 

SUBROUTINE CALCIN(IOK,ICALL,IRT,TCOM) 


CCALCIN *****SUBROUTINE CALCIN - N-LAYER ELASTIC SYSTEM **#*¥*#** 


1 


2 


ONL 


DIMENSION RR(20) ,22(20),2(5)-V (5) ,HH (4) -H (4) -AZ(190) -A(190,5) » 
TBC O FS) ee C9055) eh DiGio0,5) , Ad (190) -RJ1(190) , 800 (190) 
COMMON RR,ZZ,E,V,HH,t,yAZyAyB,CyD,Ad, RJ1, 250 
COMMON R,Z,AR,NS,N,L,ITN,P,RSZ,RST,RSR,RTR,ROM,RMU,SF 
COMMON CSZ,CST,CSR,CTR, COM, CMU 
DIMENSION TITLE(18), BZ(100), X(5,4,4%), SC(4), PM(4,4,4), FM(2,2), 
1 TEST (11) 
COMMON TITLE, PSI, NLINE, NOUTP, NTEST, TEST, 
1 ITNG, LC, JT, TZZ, PR, PA, EP, T1P, T1IM, T1, T2, T3, T4, 
2 ETO VLCC LAE ec W A pee mADUOeeD iy ak, S2l, Sad, PM, SGi, SG2, 
3 PH, PH2, VK2, VKP2, VK4, VKP4, VKK8, X, SC, FH 
DIMENSION W(4) 
W(1) = 0.34785 
w(2) = 0.652145 
w(3) = ¥(2) 
w(4) = W(1) 
VL=2.0*V(L) 
EL=(1.0+V(L)) /E(L) 
VL1=1.0-VL 
CSZ=0.0 
CST=0.0 
CSR=0.0 


Cs 


NTS1 = NTEST + 1 


ARP = AR 
IF (NOUTP) 4,4,5 


ARP = ARP*PSI 
CONTINUE 
DO 40 I=1, ITN 
INITIALIZE VHE SUB-INTEGRALS 
RSZ=0.0 
RST=0.0 
RSR=0.0 
RTR=0.0 
ROM=0.0 
RMU=0.0 
COMPUTE THE SUB-INTEGRALS 
K = 4*(I-1) 

DO 30 J=1,4 

Fi wnateas 3 
P=AZ (J) 


EP=EXP (P¥*Z) 
T1=E(J1,L) *EP 
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35 
36 
33 
40 
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T2=D(J1,L) /EP 

T1P=11+T2 

TIM=T1-T2 

T1=(A(J1,L) +B (J1,L) *Z) *EP 

T2= (C(J1,L) +D(J1,L) #2) /EP 
T2P=P*(T1+T2) 

T2M=P* (T1-T2) 

WA=AJ (J 1) *W (J) 

IF (R) 20,20,15 

BJ1=RJ1(J1) #P 

BJO=RJO (J1) *P 

RSZ=RSZ+WA*P*BJO* (VLIXTIP-T2M) 
ROM=ROM+WA*EL *BJ0* (2.0 *VL1*TIM-T2P) 
RIR=RTK+WA*P*BI1¥* (VL®T1M+T2P) 
RMU=RMU+WA*EL*BJ 1* (T1P+T2M) 
RSR=RSR+WA* (P*BJO* ( (1. 04+VL) *T1P+T2M) —BJ1* (T1P+T2M) /R) 
RST=RST+WA* (VL*¥2*BJO*¥T1P+BJ 1* (T1P+T2M) /R) 
GO TO 30 

SPECIAL ROUTINE FOR R = ZERO 

PP=P*pP 

RSZ=RSZ+WA*PP* (VL1*¥T1P-T2M) 
ROM=ROM+WA*EL*D* (2. 0*VL1*T1N-T2P) 
RST=RST+WA*PP* ((VL+0.5) *T1P+0.5*T2M) 
RSR=RST 

CONTINUE 


SP = (AZ(K+4) - AZ(Kt1))/1.7222726 
CSZ=CSZ+RSZ*SF 
CS1=CS1+RST*SP 
CSR=CSR+RSR*SF 
CTR=CIR+RTR*SF 
COM=CCM#ROM*SF 
CHU=CMU+RMU*SP 

RSZ = 2.0*RSZ*AR*SP 
TESTH = ABS (RSZ)-10.0**(-4) 
IF (I1S-NTS1 ) 31,32,32 
CONTINUE 
TEST(IIS) = TESTH 
ITS = ITS+1 
GO TO 40 
CONTINUE 
TEST(NTS1) = TESTH 
DO 33 J = VNTEST 
IF (TESTH-TEST(J)) 35,36,36 
CONTINUE 
TESTH = TEST (J) 

CONTINUE 

TEST(J) = TEST (J+1) 
CONTINUE 

IF (TESTH) 50,50,40 
CONTINUE 

IPi= i 

CS2=CSZ*ARP 
CSI=CST*ARP 
CIR=CTR*AKP 
CSR=CSR*ARP 
CCOM=CCM*ARP 
CMU=CMU*ARP 
KSTN=(CSR-V(L) *(CST#¢CSZ) ) /E(L) 
BSTS = CSZ+¢CST+CSR 
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ES 


BSTN = BSTS * (1.0-2.0*V(L))/E (L) 

ERett2a 72,7257 1 

2 = -Z 

CONTINUE 

WRITE (6,315) R,Z,CSZ,CST,CSR,CTR,BSTS,COM,RSTN, BSTN 
FORMAT (lH ,F6.2, F702, 2X, 5E13.4, 3X ,3E 13.4) 

NLINE = NLINE + 1 

IF (JT) 99,99,60 

WRITE (6,310) £ (NS) 

FORMAT (1H¢,1Z26X,4HSLOW) 

FORMAT (1+, 124X,F7.0) 

IF(Z.LE.0.00001) TCOM=TCOM+¢COM 

RETURN 

END 

SUBROUTINE PART 

#*#*e4SUBPOUTINE PART - N-LAYER ELASTIC SYSTEM ***4e4® 

DIMENSION RR(20) ,2Z(20),E(5),V (5) ,HH(4),H (4) -AZ (190) ,A(190,5), 

1 B(190,5), €(190,5) ,D(190,5) ,AJ (190) ,RJ1(190) , RIO (190) 

COMMON RR,ZZ,E,V,HH,H,AZeA,ByC,D, Ad, RI1,R250 

COMMON R,Z,AR,NS,N,L,ITN, P,RSZ,RST,RSR,RTR,ROM,RMU,SF 

COMMON CSZ,CST,CSR,CTR,COM,CMU 

DIMENSION TITLE(18), BZ(100), X(5,4,4), SC(4), PM(4,4,4), PM(2,2), 
1 TEST(11) 

COMMON TITLE, PSI, NLINE, NOUTP, NTEST, TEST, 

1 ITNG, Loo 2 lee DR Pe, EP, T1P, TIM, T1, T2, T3, T4, 
Cet oame Gwar et, 7 RA Bots BOC, Bln ih en Sale Sills PMs SG1, SG2, 
3 PH, PH2, VK2, VKP2, VK4, VKPU, VKK8, X, SC, FM 

** COMPUTE ZEROS OF J1(X) AND JO(X). SET UP GAUSS CONSTANTS ** 


BZ(1) = 0.0 
BZ(2) = 1.0 
BZ(3) = 2.4048 
BZ(4) = 3.8317 
BZ(5) = 5.5201 
BZ(6) = 7.0156 
K = ITN+1 
DO 2 I=7,K,2 
Tse LZ 


TD = 4.0*T - 1.0 

BZ(I) = 3.1415927*(T - 0.25 + 0.050661/TD 
1 -0.053041/TD**3 + 0.262051/TD**5) 

DO 3 I=8,ITN,2 

{ = (I-2)/2 

TD = 4.0*T + 1.0 

BZ(I) = 3.1415927*(T + 0.25 - 0.151982/TD 
+ 0.015399/TD**3 ~ 0.245270/TD**5) 
G1=0. 8611363 
G6 2=0. 3399810 
ZF = AR 
NIEST = 2 
IF (R) 8,8,9 
CONTINUE 
NTEST = AK/R + .0001 
IF (NIEST) 6,6,5 
CONTINUE 
NTEST = R/AR + 20001 
ZF = R 
CONTINUE 
NTEST = NTEST ¢ 1 
IF (NTEST-10) 8,8,7 
CONTINUE 
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NTEST = 10 
8 CONTINUE 
46 ** COMPUTE POINTS FOR LEGENDRE-GAUSS INTEGRATION ** 
K = 1 
ZF = 2.0*2F 
SZ2 = 0.0 
DO 28 I=1,ITN 
$Z1 = SZ2 


$Z2 = B82 (1+1)/ZF 
SReeeo 722 a SZct 
PX = SZ2 + SZ1 
SG1=SF¥*G1 
SG2=SF*G2 
AZ(K) =PX-SG1 
AZ (K+1) =PX-SG2 
AZ (K+#2) =PX*SG2 
AZ (K+3) =PX+SG1 
K=K*«4& 


28 CONTINUE 
4Q RETURN 


END 
SUBROUTINE COEE (KIN) 


CCOEE *4*4***eSUBROUTINE COEE - N-LAYER ELASTIC SYSTEM *****%* 


DIMENSION RR(20) ,2Z(20) -E(5) /V (5) ,HH (4) »H (4) -AZ(190) ,A(190,5) , 

1 B(190,5), C(190,5) -D(190,5) , Ad (190) ,RJ1(190) , RIO (190) 

COMMON RR,ZZ,5,V,Hi,HyAZ,A,B,CyD, AJ, RJ1,RI0 

COMMON R,Z,AR,NS,N,L,ITN,P,RSZ,RST,RSR,RTR,ROM,RMU,SF 

CCMMON CSZ,CST,CSR,CTR,COM, CMU 

DIMENSION TITLE(18), BZ(100), X(5e4-e4),4 SC(4), PM(4,4,4) , FM (2,2), 


1 TEST (11) 
COMMON TITLE, PSI, NLINE, NOUTP, NTEST, TEST, 
1 ITNG, POT, (T2757 0R. PA, EP, RIPsaTiIM, T1,) 22>. 63, 75; 


25 TGs 2p.) TOM, WA, BJ1, BJO, BZ, ZF, S21, SZ2, PM, SG1, SG2, 
3 PH, PHO, VK 25. UKE 2. YK4, VKPG, VKK8, X, SC, FM 
DIMENSION SV1(4,2) ,CV1(2,1) ,SV2 (4,4) »CV2 (2,2) ¢SV3 (4,8) -CV3(2,4), 
1 SV4 (4,16) ,CV4 (2,8) T(8) ,NT(4) 
COMMON SV1,CV1,SV2,CV2,SV3,CV3,SV4,CV4,T,NT 
LC = KIN 


CS-MX SET UP MATRIX X =DI*MI*KI*K*M*D 


COMPUTE THE MATRICES X (K) 
DO 10 K=1,N 

TI=E(K) *(1.04V (K#1)) /(E(K41) *(1.04V (K))) 
TIM=T1-1.0 

PH=P*H (K} 

PH2=PH*2.0 

VK2=2. 0*V (K) 

VKP2=2.0*V (K+1) 
VK4=2.0*VK2 
VKP4=2.0*VKP2 
VKK8=8. 0*V (K) *V (K+1) 


X(K,1,1)=VK4-3.0-T1 

X (K, 2, 1) =0.0 
X(K,3,1) =TIM* (PH2-VK4+1.0) 
X(K,4, 1)=-2.0*TIM*P 


T3=PH2* (VK 2-1. 0) 
T4U=VKK8+1 20-3. O*VKP2 
TS=PH2Z* (VKP2-1.0) 
T6=VKK841.0-3. OX VK2 
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X (Ky 1,2) =(T3+T4-T1* (T5476) ) /P 
X(K,4, 2) =T1M* (1. 0-PH2-VKP4) 


X (Kp 3,4) =(T3-T4-T1* (T5<T6) ) /P 


T3=EFH2*PH-VKK8+¢1.0 
T4Y=PH2* (VK2-VKP2) 


X (Ky, 1,4) = (T3+T4+VKP2-T1* (T3+T4+VK2)) /P 
X(K,3e2)=  (~T3+T4-VKP2+T1* (T3-T4+4VK2)) /P 


X (Kg 1,3) =T1M* (1. O-PH2-VK4) 
X (K, 2,3) =2-0*T1M*P 

X (K,3,3) =VK4-3.0-T1 
X(K,4,3) =0.0 


X (K, 2,4) =T1M* (PH2-VKP 441.0) 
X(K,4,4) =T1* (VKP4-3.0) =1.0 
K = K 

CONTINUE 

COMPUTE THE PRODUCT MATRICES PM 
SC (N) =4.0* (V(N)-1.0) 

IF (N-2) 13,11,11 

DO 12 K1=2,N 

M=NS-K1 

SC (M) =SC (M+#1) *4. 0* (V(M) -1.20) 
CONTINUE 

CONTINUE 


DO 14 J=1,2 

S Wal bs was (Ky 1, 0+ 2) 
CONTINUE 

CV1(1,1) = -2.0*P*H(K) 
CV1(2,1) = 0.0 

K = K-1 

IF(K) 50,50,20 


CONTINUE 
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SV1(4,J) 
DO 21 I=1,4 
SV2(I,J1-1) = X(K,I,1) *T(1) +X (K,1I, 2) *T (2) 


SV2(I,J1) = X(K,I,3) *T(3) +X (K,I,4) *T(4) 
CONTINUE 

TEV) = GVA(1, 1) 

T(2) = -2.0*P*H(K) 

CV2.( tat (1) 

CV2Z(1,2) = T(2) 

CV2.(2¢1) = 1(1)-7 (2) 

CV2(2,2) = 0.0 

K = K-1 


IF (K) 50,50,30 
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CCNTINUE 

DO 34 J=1,4 

TUG 2 

DeItat-2) 32,32,31 

Ji = 3142 

CONTINUE 

T (1) SV2(1,J) 

T (2) SV2 (2,J) 

T (3) $V2(3,J) 

T (4) SV2 (4,J) 

DO 33 I=1,4 

S¥3 (I,d1) = X(KoI_1) *T (1) +X (KeIe 2) *T (2) 
SV3(I,J1+2) = X(K,I,3) *T(3) +X (K,1,4) *T(4) 
CONTINUE 

T(1) = -2.0*P*H(K) 

DO 35 J=1,2 

CV3 (1,3) = CV2(1,J) 

CV3 (2,3) = CV2 (1,3) -T (1) 

Cvsemu+2) = CV2 (2,5) +71) 


tf nw tou 


CV3(2,J3+2) = CV2(2,J) 
CONTINUE 

K = K-1 

IF (K) 50,50,40 
CONTINUE 

DO 42 J=1,4 

T(1) = Sv3(1,J) 
T(2) = SV3(2,d) 
T(3) = SV3(3,d) 
T(4) = SV3(4,d) 
T(5) = SV3(1,Jd+4) 
7 (6) = SV3(2,J+4) 
T(7) = SV3(3,J3+4) 
T(8) = SV3(4,J+4) 


DO 41 I=1,4 

SVG (I,I) = X(KyI, 1) *T(1) +X (K,I,2) *2(2) 
SV4 (1,d+4) X (K,1I,3) *T (3) +X (K,1,4) *T (4) 
SV4(I,d+8) = X(K,1I,1)*T(5) ¢X(K,I, 2) *T (6) 
SV4(I,d+12) = X(K,I,3) *T(7) +X (K,I,4) *T (8) 
CONTINUE 

T(1) = -2.0*P*H(K) 

DO 43 J=1,4 

Cv4 (1,0) 7 S09 CV 31,3) 

CV4 (2,5) = CV3(1,J3)-T(1) 

CV4 (1,9 +4) CV3(2,J) +T(1) 

CV4(2,0+4) = CV3 (2,4) 

CONTINUE 
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CONTINUE 

NT(1) = 1 

DO 51 K=2,N 

NI(K) = NT (K-1)+NT (K-1) 
DO €0 K=1,N 


K1 = NS-K 

DO 52 I=1,4 
PM(K1,1I1,1) = 0.0 
PM(K1,I,2) = 0.0 
CONTINUE 


I1 = NT (ik) 
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DO 80 I=1,11 
I2 = I+1I1 


GO 10 (61,62, 


CONTINUE 
T (3) 


GO 10 65 
CONTINUE 


Taye = cv 2 C1:, 
Tie = CV 2 (2, 


GO TO 65 
CONTINUE 
T (3) 


GO 10 65 
CONTINUE 


T(3) = cV4(1, 


T(4) = CV4(2, 
CCNIINUE 

T(1) = 0.0 
T(2)-= 0.0 


= cv1(1, 
T(4) = CV1(2, 


= CV3 (1, 
T(4) = CV3(2, 
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IF (T(3)+68.0) 67,66,66 
T(1) = EXP (T(3)) 
IF (T(4)+68.0) 69,68,68 
T(2) = EXP (T(4)) 


CCNIINUE 

DO €0 J=1,2 

GO! TO) (31072573, 79)_.K 
CONTINUE 

T(3) = SV1(J,I) 
T(4) = SV1(J,12) 
@(5) = SV1(J*2,1) 
T(6) = SV1(J3+2,12) 
GO T0 75 

T(3) = SV2(J,1) 
T(4) = SV2(J,12) 
T(5) = SV2(J+2,1) 
T(6) = SV¥2(J+2,12) 
GO TO 75 

T(3) = SV3(J,I1) 
T(4) = SV3(J,I2) 
T(5) = SV3(J+2,1) 
T(6) = SV3(J+2,1I2) 
GO 10 75 

T(3) = SV4(J,I) 
T(4) = SV4(J,12) 
T(5) = SV4(J+2,1) 
T(6) = SV4(J+2,12) 
CONTINUE 


PM(K1,J,1) 
P&(K1,J,2) 
PM(K1,J+2,1) 
PM(K1,J+2, 2) 
CCNTINUE 


PM(K1,J, 1) #2(1) *T (3) 
PM(K1,J,2) #1 (1) *T (4) 
= PM(K1,d+2,1) +T (2) *T (5) 
= PM(K1,J+2, 2) +T (2) *T (6) 


SCLVE FOR C(NS) AND D(NS) 


V2=2.0*V (1) 
Vz1=Vz-1.0 
DO 90 J=1,2 


PM(1,J)=P*PM(1,1,0) #V2*PM(1,2,d) ¢P¥PM(1,3,5)-V2*PM(1,4,d) 
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90 FM(2,J) =P*PM(1,1,J) +¥21*PM(1,2,J)-P*PM(1,3,d) #+V21*PM(1,4,J) 
DPAC=SC (1) / ((FM (1,1) *FM (2,2) -FM(2, 1) *FM(1,2)) *P*P) 


A(LC,NS) = 0.0 
B(LC,NS) = 0.0 
C(LC,NS) = -FM(1,2) *DFAC 
D(LC,NS) = FM(1,1)*DFAC 


C BACKSOLVE FOR THE OTHER A,3,C,D 
DO 91 K1=1,N 
A(LC,K1)=(2“(K1, 1, 1) *C (LC, NS) #PM(K1,1,2) *D(LC,NS) ) /SC (K1) 
B(LC, K1)= (2M (K1,2,1) *C(LC, NS) ¢PM (K1, 2,2) *D(LC,NS) ) /SC (K1) 
C (LC, K1) = (24 (K1, 3,1) *C (LC,NS) ¢PM(K 1,3, 2) *D (LC, NS)) /SC (K1) 

91 D(LC,K1)=(PM(K1, 4,1) *C(LC,NS) ¢PM (K1,4,2) *D (LC,NS) ) /SC (K1) 

RETURN 
END 
SUBROUTINE BESSEL(NI, XI, Y) 

CBESSEL #***##SUBROUTINE BESSEL - N-LAYER ELASTIC SYSTEM ###¥4x 
DIMENSION PZ(6) »QZ(6) »P1(6) -Q1(6) ,D(20) 


c 
Cc 
- 1 PZ(1)=1.0 
PZ(2) = -1.125E-4 
PZ(3) = 2.87109E-7 
PZ(4) = -2.34495-9 
PZ(5) = 3.98068E-11 
PZ(6) = -1.15361E-12 
c 
QZ(1) = -5.0E=3 
QZ(2) = 4.6875E-6 
QZ(3) = +2.32559E-8 
QZ(4) = 2.83071E-10 
QZ(5) = -6.39121E-12 
QZ(6) = 2.31247z-13 
S 
Cc 
aw by aa 
P1(2) = 1.87S5E-4 
P1(3) = -3.69141E-7 
P1(4) = 2.77132E-9 
P1(5) = -4.51144E-11 
P1(6) = 1.27505E-12 
Cc 
Q1(1) = 1.5E<2 
Q1(2) = -6.5625E-6 
Q1(3) = 2.842385-8 
Q1(4) = +3, 266205E-10 
Q1(5) = 7.14312E-12 
Q1(6) = -2.53271E-13 
C 
Cc 
PI = 3.14159 
PI2 = 2.0*PI 
Cc 
- 
9 : = NI 
= XI 
_ (X-7.0) 10,10, 160 
Cc 
10 X2=%72.0 
FAC=-X2*X2 
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T=I 
C=FAC*Cy (T*T) 


TEST=ABS (C) — 10.0** (-8) 
LE CEES) mahie gt idee 


Y=Y+C 
CONTINUE 
C=X2 

Y=C 

DO 16 I=1,34 
T=I 


C=FAC*C/ (T* (T+#1.0)) 
TEST=ABS (C) - 10.0**(-8) 
IF (TEST) 17,17,15 


Y=Y+tcC 
CCNTINUE 
RETURN 


IF (N) 161,161,164 


DO 162 I=1,6 
D(I) = PZ(I) 
D(I+10) = Q2Z(I) 
CONTINUE 

GO 10 163 


DO 165 I=1,6 
D(I) = P1(I) 
D(I+10) = Q1(I) 
CCNIINUE 
CONTINUE 

T1 = 25.0/X 
T2=T14T1 

P = D(6) *T2+D(5) 
DO 170 I=1,4 


Je= 5-1 
P = P*T2+D (J) 
CONTINUE 


Q = D(16) *T2+D (15) 


DO 171 I=1,4 


J = 5-I 

C = Q¥*T2+D(J+10) 
CONTINUE 

Q = Q*Ti 

T4 = SQRT (X*PI) 
T6) =<SIN (X) 

T7: =ecOSs (x) 


IF (N) 180,180,185 


LS = iC (P=0) *T6 + 
GO TO 99 
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